NPS  ARCHIVE 

1961 

LINDER,  I.                         [i|| 

k)u!hl 


. 


ill 


(it 


COR 

R 

li 

JBHHI 

uiK/1  i;    i 


•  ■ ■  >.\' 1 1 


ELATION  PROCESSES  IN 
ANTENNA  ARRAYS 


it  (Hf)iJllM  IP5  IT  r! I'fi'il 
lii'Hplilfiin  ri^lHillilinltHilHiufiiiillHK 

iir In P(l!$'  t  IH I iiH,iH'J'n  i!r 
|l llillill'iili! IwiikWili  ilm'l  i  il^' .  1 1 
JUIhU  luwuhlflimo  ml  111  III  ll  ImilluN 


Hi!         1*   ?t 

lifl  M   lU    I    I 


it  Qui 


lipi 


I  liiHiUf 


mill 
jfillllim  It 


I 
1 

| 

Jl«  Ml 

Mil 

m 


I 

1 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOt 
MONTPflEY  CA   93943-5101 


/?6/ 


Correlation  Processes  in  Antenna  Arrays 

By 

Isham  Wiseman  Linder 
// 

3.S.  (United  States  Naval  Acadeny)  I9I46 

B.S.  (Massachusetts  Institute  of  Technology)  19U9 

M.S.  (United  States  Naval  Postgraduate  School)  1956 

DISSERTATION 

Submitted  in  partial  satisfaction  of  the  requirements  for  the  degree  of 

DOCTOR  OF  PHILOSOPHY 

in 
Engineering  Science 

in  the 
GRADUATE  DIVISION 
of  the 
UNIVERSITY  OF  CALIFORNIA 


Library 

U.  S.  Navnl  Postgraduate  School 

Monterey,  California 


RRELATION  PRC  TENNA  ARRAYS 

TABLE  .      TENTS 

Introduction  -, 

II.   Antenna  Patterns  of  Nonlinear  Arrays  l 

A.   General 


B.  'chromatic  Signals  5 

C.  Band-Limited  Random  Signals ,  General  Case  10 

D.  Sand-Limited  Random  Signals,  Specific  Cases 

III    Effect  of  Finite  Averaging  Time  25 

Resolution  Characteristics  of  Correlation  Arrays  35 

Effect  of  Distributed  Noise  on  the  Nonlinear  Antenna 

stem  43 


T 


V. 


References 


52 

56 
56 


VI.   Antenna  Arrays  and  Channel  Capacity 

VII.   Design  of  Nonlinear  Arrays 

A.  Basic  Considerations 

B.  Products  of  Linear  Arrays  57 

C.  Data  Processing  by  Correlation  Circuits  59 

>pendix  I.   The  Probability  Distribution  for  the  Output 
koise  Voltage  of  a  Two- Element  Correlation 
Array  67 

Appendix  II.   Calculation  of  Voltage  Response  Patterns 

with  Two  Signal  Sources  74 

Appendix  III.   Application  of  Correlation  Techniques  to 

Linear  Arrays  01 


108 


CORRELATION  PROCESSES  IN  ANTENNA  ARRAYS 

I.   INTRODUCTK 

The  design  of  an  antenna  array  whose  performance  must 
be  optimized  in  some  specific  sense  (maximum  signal  to  noise 
ratio  in  the  output  circuit,  minimum  antenna  pattern  beamwidth, 
etc.)  when  the  array  is  to  receive  a  signal  buried  in  a  general 
noisy  medium  can  be  accomplished  through  noise  theory  con- 
siderations and  correlation  techniques.   These  statistical 
methods  can  be  applied  to  two  general  types  of  multi- element 
antenna  arrays.   First,  there  is  the  normal  linear  array  in 
which  the  output  circuit  contains  the  sum  of  the  voltages  in- 
duced on  the  individual  elements.   Analysis  of  the  noise  volt- 
ages in  terms  of  their  correlation  coefficients  will  lead  to 
a  specification  of  element  spacing  which  will  optimize  the 
antenna  performance.   Second,  there  is  the  newer  nonlinear 
array  in  which  both  multiplication  and  summation  processes 
may  be  used  in  combining  the  individual  element  voltages  to 
form  a  single  array  output  voltage.   This  type  of  array  re- 
sults from  an  attempt  to  implement  the  definition  of  the  cor- 
relation coefficient  by  physical  circuitry,  and  can  be  ex- 
tensively analyzed  by  methods  of  statistical  noise  theory. 

Correlation  techniques  have  been  introduced  into  the 
analysis  of  the  linear  array  primarily  by  acoustical  engineers" 


1,2References  are  listed  at  the  end  of  this  report. 


2. 

The  manner  in  which  the  correlation  coefficient  enters  into  a 
consideration  of  this  general  type  of  array  can  be  seen  from 
the  simple  example  of  the  two- element  linear  array.   A  random 
noise  source  will  induce  a  voltage  on  one  element  which  will 
be  identical  to  that  induced  on  the  other  element  except  for 
a  time  delay  caused  by  the  time  of  travel  of  the  wavefront 
between  elements.   If  e-^(t)  and  e2(t)  represent  the  instan- 
taneous voltages  on  the  two  elements,  then  the  average  power 
in  the  output  circuit  of  the  array  becomes: 

ez  =  [ei(t1)  +  e2(t2)]z 


=  ei(tL)2  +  e2(t2)2   +  2  e1(t1)e2(t2) 
*  e1(t1)z  ♦  e2(t2)z  +  2  R12(t1-t2) 

where  R-.2(^f)  is  the  cross-correlation  coefficient. 

The  cross-correlation  coefficient  will,  in  general,  take 
on  both  positive  and  negative  values  as  «  is  varied.   Thus 
it  should  be  possible  to  select  an  element  spacing  which  will 
minimize  the  average  noise  power,  e- .   This  rather  simple  idea 
can  be  applied  to  multi-element  linear  arrays  --  the  output 
voltage  is  expressed  as  the  sum  of  noise  power  contributions 
of  each  element  plus  the  cross  correlation  coefficients  of  all 
pairs  of  elements  in  the  array.   This  total  average  noise 
power  is  then  minimized  by  proper  spacing  of  the  antenna 
elements. 

Some  previous  results  in  this  application  of  correlation 
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techniques  to  linear  arrays  are  presented  in  Appendix  I,  to- 
gether with  a  general  methodology  for  applying  these  principles 
to  the  types  of  antennas  useful  in  electromagnetic  radiation. 

The  main  body  of  this  report  will  be  concerned  with 
the  analysis  of  the  second  type  of  array- -the  array  which 
utilizes  multiplicative  processes  in  combining  the  individual 
element  voltages  to  form  the  single  array  output  voltage.   The 
initial  section  will  be  concerned  with  a  review  of  antenna 
response  patterns  which  result  from  these  nonlinear  processes 
in  the  antenna  circuitry.   Succeeding  sections  will  discuss 
the  effects  of  distributed  noise,  of  multiple  signal  sources, 
and  of  a  variable  averaging  time  on  the  output  voltage.   Finally 
the  results  of  these  analyses  will  be  utilized  to  reach  general 
conclusions  about  the  type  of  nonlinear  antenna  which  will  give 
optimum  performance,  and  examples  of  these  antennas  will  be 
developed. 


4. 

II.   ANTENNA  PATTERNS  OF  NONLINEAR  ARRAYS 
A.  General 

The  principle  of  multiplying  together  the  output  volt- 
ages from  two  antennas  is  well  known;  its  applications  in- 

3  4  5 
elude,  for  example,  radio  astronomy.   Several  papers  J*^*J   have 

discussed  some  aspects  of  the  subject,  but  no  treatment  appears 
to  have  been  made  of  the  problem  of  development  of  nonlinear 
arrays  intended  to  detect  a  signal  buried  in  a  medium  of  ran- 
domly varying  noise  sources  or  a  signal  partially  hidden  by 
other  undesired  signals. 

Before  discussing  the  types  of  antenna  directional  pat- 
terns resulting  from  the  use  of  nonlinear  processes  in  the 
antenna  circuitry,  it  would  be  best  to  define  the  general  type 
of  antenna  array  which  will  be  employed  and  to  establish  the 
general  description  of  the  randomly  varying  noise  and  signal 
sources. 

To  remain  completely  general,  it  will  be  assumed  that 
the  noise  and  signal  voltages  will  be  applied  to  the  M  elements 
of  a  receiving  array,  and  that  these  elements  are  connected 
to  individual  amplifiers  (A^,  A2 ,  -  -  ->am^  an<*  ph^se  shifters 
(01  >  $2  >  "  "  "»  ^M^  *   ^e  outPut  °f  each  antenna-amplifier- 
phase  shifter  circuit  is  connected  to  a  combining  network 
which  will  operate  on  the  individual  element  voltages  to  pro- 
vide an  optimum  output  voltage  from  the  total  array. 
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The  combining  network  will  consist  of  various  summation 
and  multiplication  circuits.   For  practical  effectiveness  time 
averaging  will  be  required  in  this  network.   The  combination 
of  multiplication  and  time  averaging  then  satisfies  the  defi- 
nition of  the  correlation  function: 

T 
R..(t)  =  litn  1  \  n.(t)  n.(t  + 1)  dt. 

In  considering  practical  networks  that  can  be  used  to 
carry  out  this  correlation  process,  two  general  variables  arise 
(1)  possible  combinations  of  the  multiplication  and  averaging 
circuits,  and  (2)  restriction  of  time  averaging  to  a  finite 
time  interval.   Both  of  these  will  be  considered  in  subsequent 
sections,  and  will  form  the  basis  for  the  final  statement  of 
the  optimum  type  of  nonlinear  array. 

For  this  analysis  assume  that  the  noise  sources  are 
spherically  distributed  and  are  statistically  independent. 
The  voltage  produced  on  an  isotropic  antenna  by  the  noise 

originating  in  an  element  of  solid  angle  dXJL  can  be  repre- 

f\    7 
sented  by  the  usual  Fourier  series:  ' 

OP 

x(t  )   *^jT     [ak   cos(2-rr^t)  +  bksin(2nfkt)      ] 
k=o 


6. 


where 

T^tD        T'Ei  ^^1    =  °'  for  m  #  n 

lira  T-E{_bmbn\    -  0,  for  oi  *  n 


lira  f'Eia  b  I 


lira  T.E{an     }  -ll«        T-E^2]    -   2  W(fn) 


2 


In  this  representation,  the  coefficients  a.  and  b,  are  dis- 
tributed normally  with  means  zero;  E  is  the  expectation 
operator,  and  W(f)  is  the  (two-sided)  power  spectral  density 
of  the  voltage  induced  on  the  antenna  elements  by  the  source. 
It  should  be  noted  that  the  inherent  frequency  response  of  the 
antenna  elements  and  amplifier  circuits  will  lead  to  the 
familiar  narrow-band  gaussian  random  process. 

The  directivity  pattern  of  the  antenna  array  will  be 
defined  to  be  the  vol  Cage  response  at  the  output  of  the  an- 
tenna array  for  an  arbitrary  angular  location,  0,  of  a  single 
signal  source;  for  simplicity  this  will  be  the  two-dimensional 
pattern. 


XT     V     XT  17    t? 

Line    of    At«*x 

B.   Monochromatic  Signals. 

The  antenna  directivity  patterns  which  result  when 
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individual  element  voltages  induced  by  a  single  monochromatic 

source  (negligible  noise)  are  correlated  have  been  investigated 

3  4 
by  several  researchers.  »    In  general,  the  procedure  has  been 

to  recall  the  directivity  patterns  of  linear  additive  arrays 
and  then  to  demonstrate  product  arrays  which  will  give  equiva- 
lent directivity  patterns. 

The  directivity  patterns  for  linear  arrays  with  constant 
element  spacing,  d,  can  be  expressed  by  one  of  two  series,  de- 
pending on  an  odd  or  an  even  number  of  elements  in  the  array. 

1.  For  an  array  of  2n+l  elements,  the  sum  pattern  is 

given  by 

n 

P2n+l(t>  =  £_    \  cos  2ku 
k=0 

where  u  =  (ird/jL  )  sin  4>.      The  array  is  assumed  to  be  sym- 
metrical with  respect  to  its  center  element,  A_,  (A.  =  A  .). 

0    i    -i 

P2n+^(^)  is  the  amplitude  gain  given  by  the  array  to  an  input 
signal  Eq  sin  (<Ot+VO.  The  center  of  the  array  is  the  phase 
reference. 

2.  A  linear  array  of  2n  elements  and  constant  spacing, 
d,  has  a  sum  pattern 

n 
P2nW  =^T  Bk  cos(2k  -  l)u 

k=l 
Again  symmetry  is  assumed,  so  3^  is  the  common  amplification 

of  the  two  elements  which  have  distance  (k  -  £)  d  from  the 

2 
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array  center. 

Berman  and  Clay  have  described  an  array  in  which  pairs 
of  signals  are  multiplied  together  and  then  time  averaged. 
The  output  voltage  is  formed  by  multiplying  together  all  of 
the  resulting  time  averages.   With  this  procedure,  an  array 
of  n  +  1  elements,  with  successive  elements  spaced  D,  2D,  4D, 
.  .  .  ,  2n"  D  from  the  first  element,  has  a  directivity  pattern 
equivalent  to  that  of  a  2n  element  additive  array  with  constant 
element  spacing  2D.   In  this  case  the  product  array  is  about 
one- fourth  the  length  of  the  equivalent  linear  array.   Such  a 
product  array  with  four  elements  could  be  represented  as 

\ 


»--f »[x] >  AV<S* 

—4 J  [* 

o-~) >   X  I »|  AVG.| ' 


3* 

4< 


*  x 


*  AVG, 


l 


X 


4  fT 

?T"   lim  21     J 
1     !  T-*oo     "'-T 


V    eel    lim 

i=2 


cos*)t  co  s  ( co  t-«*t  ^  )  dt 


V    oc 
out 


n 

i=2 


cos 


*>T 


li 


With  proper  spacing  of  the  elements,  this  pattern  is  the  same 
as  that  of  an  eight-element  linear  array. 
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A  second  possible  product  array  carries  out  all  desired 
multiplications  before  finally  time  averaging  the  product. 
Analysis  of  the  directivity  pattern  in  this  case  shows  a  smaller 
saving  in  overall  length  over  the  equivalent  linear  array  than 
was  the  case  in  the  first  product  array  discussed.   Again  a 
four- element  example  is 


I  o- 


2* 


3- 


40. 


^ut00  ltm 
T-oco 


r 

-T   i=l 


cos(o*  t  ~o/t  .  )dt 


With  proper  spacing  of  the  elements,  this  pattern  is  the  same 

as  that  of  a  six- element  linear  array. 

It  is  also  possible  to  start  with  two  elements  and  to 

process  the  element  voltages  to  get  directional  patterns  which 

are  equivalent  to  those  of  odd  or  even  element  additive  arrays, 

i  # 


i 


z 
o 


x  k- 


AVQ.  I 


10. 

Voltage  induced  on  element  #1     sin(tt# t  +^) 

Voltage  induced  on  element  #2     sin[«Mt  -  (*£  sin^O  +|M 

Output  of  multiplier       Y=  sin(wt  +yOsittiJMt~  (££  sin  $0  +*] 

After  averaging  T"=  2  cos  &*     sin  ^ 

-   1     ird  1 

And,  letting  d  =  2s,  Y  =  -  cos  (if  sin  4  )  K   -  cos  u 

Brown  and  Rowlands  have  shown  how  to  use  this  function 
to  synthesize  a  directivity  pattern.   Using  the  relation 

T  (cos  u)  =  cos  nu         n=0,l,2,  ... 

where  T  (x)  is  the  Chebyshev  polynomial  of  the  first  kind,  every 
term  in  the  cosine  series  of  the  desired  directivity  pattern 
can  be  synthesized  by  a  suitable  combination  of  power  law  de- 
vices operating  on  the  output  of  the  two- element  correlation 
pair.   Arithmetically,  two  elements  can  thus  be  made  to  give 
a  directivity  pattern  equivalent  to  that  of  a  linear  additive 
array  of  any  arbitrary  length.   In  practice,  noise  consider- 
ations, ignored  in  these  pattern  calculations,  and  the  presence 
of  more  than  one  signal  source  would  limit  the  length  of  the 
equivalent  linear  array  which  could  be  synthesized.   These 
points  will  be  considered  in  more  detail  in  a  later  section. 

C.   Band-limited  Random  Signals,  General  Case 

As  with  the  monochromatic  signals,  the  basic  unit  in 
the  study  of  correlation  of  band- limited  random  signals  is  the 
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two-eleraent  array.   In  this  case  the  signal  can  be  described 
only  in  statistical  terms,  and  statistical  methods  can  be 
utilized  to  describe  the  performance  of  the  system.   As  be- 
fore, noise  will  be  neglected  in  initial  considerations  and 
will  enter  in  a  later  section.   Correlation  in  such  a  two- 
eleraent  array  is  directly  analogous  to  the  operation  of  cor- 
relation detectors,  which  have  been  described  in  the  litera- 
ti 
ture.   Performance  of  correlation  detectors,  however,  is 

usually  studied  in  terms  not  directly  applicable  to  the  des- 
cription of  the  characteristics  of  antenna  arrays. 


lo 


*,<tV 


V.tr) 


2* 


XiCt). 


9ctT 
6«^ 


Vf<t) 
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Y(t) 


*  AV6. 


Yft) 
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We  assume  that  the  voltages  induced  on  the  antenna  elements 
are  stationary  and  ergodic  random  functions  of  time  with  nor- 
mal distributions  and  are  confined  to  Iff  B,  a  closed  interval 
in «a.   The  amplifiers  are  linear,  time- invariant  networks  with 
system  functions  H(u> )  =  (h(u>)|   e^*    and  G(*>  )   =|G(u>)\e* 
Fourier  transformable  into  the  impulse  responses  h(t)  and  g(t). 

For  each  set  of  the  ensemble  of  signal  functions,  the 
Fourier  expansion  is: 

M 
xx(t)  =  ^T  a1  cos(u*tt  +  ^i) 
i=l 


12. 
N 
x2(t)   ■  JT      b.   cos(4*Jit   +Jf».  ) 

where 

— <r 

1.  a.,  b.  are  Rayleigh  distributed  with  a.   *  2W.(f.)  5f 
and  b^  ■  2W2(f.)Sf,  with  W(f.)  as  the  power  density  over  a 
frequency  interval  if  centered  at  f  ; 

2.  ^-,V,  have  uniform  distributions  over  (0,  2w); 

3.  a*  is  the  lower  edge  of  the  band  of  width  B  cps,  and 
<*>L  =  *>0  +  2«i/T;  i  =  1,  ..  .,  N. 

Then 

N 


^(t)   =     £~      ai   hi   COS<u>   *  +  4^   +1^ 


i=l 

N 


y2U>  =  ]T[  b.  g.  cosc^^  +  y.  ♦  r: 


i=l 
After  multiplication 


N   N 


Y(t)  =  yi(t)y2(t)  .£  21  ^VAtH^V^Y'^V 


i=»l  j=l 


Ico.^^  -"f+4k^    %*A* 


The  first  cosine  term  describes  a  band  of  frequencies  centered 
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at  twice  the  mid- frequency  of  the  band-limited  signal,  and  the 
second  cosine  term  describes  the  dc  and  low  frequency  component 
of  the  output  voltage. 
After  averaging: 

00 
YTt7  ■  \   a  b  g  h  cos(*  -*+%    -f  )d*>. 

0 

Now  in  this  antenna  case  X2(t)  «  x  (t)  e     .    So 

a,  =  b.  and  y».  ■  4.    +  «••  H  .   And  if  the  amplifiers  are  assumed 
i    i      i    l    i  ^ 

to  have  identical  phase  functions  and  constant  amplitude  func- 
tions, the  output  is 

oo 
YTt7  =(GH)  \  W(<*>)  cos*>T  dm   =  (G  H)  R(T). 

0 
of  the  correlation  function.   For  this  normally  distributed 
signal  with  a  rectangular  frequency  function  (center  frequency 


V 


YTtT     =      (G-H)R   Sl^BY       cos   2trfr 

WBt  0 


The  obvious  advantage  to  specifying  amplifiers  with 
identical  phase  functions  and  constant  amplitude  functions  is 
that  this  leads  directly  to  a  statement  of  the  output  of  the 
system  in  terms  of  the  correlation  function  of  the  input  band- 
limited  signal.   A  signal  which  is  normally  distributed  over 
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the  pass  band  of  the  amplifiers  can  be  analyzed  most  simply 
through  its  characteristic  function.   As  a  result  of  these  two 
statements,  the  performance  of  correlation  arrays  with  more 
elements  can  be  determined  in  a  direct  manner. 

Before  considering  more  elements,  however,  it  is  useful 
to  point  out  the  self  noise  or  fluctuation  component  of  the 
output  voltage  of  the  multiplier  in  the  simple  two-element  cor- 
relator.  For  a  normally  distributed  signal  source  located  in 
the  principal  lobe  of  the  two-element  array  (T*  ■  0),  the  out- 
put of  the  multiplier  has  a  Chi-square  distribution  whose  mean 
value  is  the  desired  array  output  voltage,  and  whose  ac  com- 
ponent is  the  self  noise  which  must  be  minimized  by  time 
averaging.   For  an  arbitrarily  located  source  it  was  shown  pre- 
viously that  the  mean  value  is  proportional  to  R(T),  the  cor- 
relation coefficient.   The  variance  is  proportional  to 


[R(0)2  +  R(T)2],  or 


•2  ccR(T)2 


j"  L   +  RCO^  ] 
L      R(T)2  J 


This  variance  does  not  go  to  zero.  Its  value  decreases  by  1/2 
as  the  signal  moves  from  the  main  lobe  to  a  position  where  its 
mean  value  is  zero. 

As  an  example  of  the  calculations  of  directivity  pat- 
terns of  correlation  arrays  with  a  greater  number  of  elements, 
consider  a  four-element  array. 
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3  A 

I — m* — 1 


The  desired  output  voltage  is  simply  the  expected  value  of  the 
product  of  the  four  element  voltages. 


And,  since  these  are  normally  distributed  band- limited  voltages, 
this  becomes 

Y(t)  =^y1(t)v2(t^v3(t)v4(t^  +fv1(t)v3(tj[v2(t)v4(t^  ♦ 

[v1(t)v4(t^[v2(t)v3(t)]. 

ecR(T1?)  R(*t,,)  ♦  R(r  )  R(T  )  +  R(T  )  R(T   ). 
12      34        13      24        14      23 

With  a  rectangular  frequency  function,  each  correlation  co- 
efficient has  the  form 


R(Tij)  ■  R(0) 


sintr  BT 


ij cos  2lff_T 


nBT 


O^ij 


^j 


And  with  relative  element  spacing  chosen  to  make  &.~   m   D, 

d«o  -  D,  and  d~,  =  2D,  the  mean  output  voltage  can  be  put  in  a 
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form  resembling  the  directivity  pattern  of  a  uniform  linear 
array. 

y(t')   rr  -  lsin  BX  sin  2BX        sin   2BX  sin   3Bxl  .  foe  vN 

Y(t)  "  2  [—&        wT  +  "^ix~~   3bTJ  cos  1(2f0x) 

1  \  sin  BX  sin  2BX  .  sin  BX  sin  4Bxl.,_  ,/,-  v* 

2  [— BX 7HT-  *  ~" BX 5EX Jcos  3(2f0x)  * 

1  fsin  BX  sin  4BX   sin  2BX  sin  3BX|      -f9#  r» 

2  BX     4BX       2BX     3BX  J  COS  DU£0X; 

where  X  =  (tt  D/c)  sin  fl(.   For  a  very  narrow  bandwidth  the  co- 
efficients of  the  cosine  terms  approach  unity  and  the  direc- 
tivity pattern  of  this  correlation  array  (length  4D)  is  the 
same  as  that  of  a  six-element  uniform  linear  array  (length  10D) . 
As  the  bandwidth  increases,  these  coefficients  influence  the 
amplitude  of  the  sidelobe  structure  of  the  directivity  pat- 
tern.  This  is  shown  in  Fig.  1  for  one  basic  spacing,  D,  and 
for  two  bandwidths. 

The  assumption  of  a  random  signal  source  with  a  normal 
distribution  permits  a  pattern  analysis  by  the  characteristic 
function  method.   The  characteristic  function  is 

vr>4 j?1  Xl+j<2X2  + +if^f  -^ImtjmJ 

and  the  pattern  will  be  given  by  the  appropriate  coefficient  in 
the  expansion  of  this  function.  For  example,  consider  a  simple 
three-element  array.   The  directivity  pattern  will  be  given  by 
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iv2Vj> 


T  T  5 

I — <5> 1        .ra» 1 

2rj 


where  E  Jxjis   the  expectation  operator.      This  will  be  the  co- 
efficient of    [   (jf  l)A    .    (jf  2)2/2.      .    (jf3)/ll   ln  the   ^ 

ies   expansion  of  exp[-(l/2)iiRCr      )f     #    1 

mix  rs    *  r  "*  s J 

E  [v^fc]  oc    R(0)R(T13)   ♦  2R(T12)R<*„) 


ser 


If  we   assume   a  rectangular   frequency   function  and  relative 
element   spacings  d^   =  D,    d^   =   3D,    the  directivity  pattern 
can  be  written 


V 


cut 


-ft 


anK sin  3BX 


^cD.2(a0X)^^«j°K   4n3B^os4(2f< 


^j:os4(2f0X)j 


where  X  «  (tr  d/c)  sin  4. 

The  variance  of  the  output  voltage  is  a  measure  of  the 
unwanted  ac  power  which  must  be  reduced  by  time  averaging,  so 
this  quantity  must  be  known  to  complete  the  description  of  a 
correlation  array.   Just  as  in  the  calculation  of  the  average 
output  voltage,  the  variance  can  be  determined  by  the  charac- 
teristic function  method,  in  this  case  by  first  finding 
E  \(V1V22V3)2J  .   This  is  the  coefficient  of 
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1.  0 


0.  5.„ 


--    0.  5 


FIG.    1  - 


-  Nor malted  directivity  pattern,    four-element 
correlation  array. 
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[  (jf  t)2/2.  *  (jf  2)4/4.  ■  (jf3)2/2.]  in  the  same  power 
series  expansion. 

a|CV1V22V3)2joc  [3   R(0)4  +   6   R(0)2R(TL3)2   ♦   12R(0)2R(T  ^   ♦ 

12R(0)2R(T23)2+   24R(T12^R<r23)2+  48R(0)RCr12)R<r2^R(r13^. 

The  mean  and  the  variance  of  this  array  voltage,  each  normalized 
to  unity  at  (2f  X)  =  0° ,  are  shown  in  Fig.  2. 


D.  Band-Limited  Random  Signals,  Special  Cases 

In  the  previous  case  no  frequency  restriction  was  im- 
posed on  the  multiplication  circuitry.   It  is  apparent  that  as 
the  number  of  antenna  elements  increases  this  operation  of  mul- 
tiplication will  impose  frequency  bandwidth  requirements  which 
may  be  difficult  to  meet.   For  each  stage  of  multiplication 
added  to  the  array  circuitry,  frequency  requirements  are 
doubled. 

As  shown  in  the  equations  for  the  two  element  product 
array,  the  output  voltage  consists  of  a  low  frequency  component 
and  a  component  centered  at  2f_.  If  this  low  frequency  com- 
ponent of  each  product  pair  is  selected  by  a  low  pass  filter 
for  use  in  subsequent  stages  of  multiplication,  the  over-all 
frequency  bandwidth  requirements  of  the  correlation  circuits  are 
kept  within  reasonable  bounds. 

In  this  case  the  directivity  pattern  cannot  be  obtained 
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fTom  the  characteristic  function  but  must  be  calculated  from 
the  low  frequency  term  of  the  equation  for  the  two-elemen  pro- 
duct pair.   By  rejecting  the  double  frequency  component  some 
increase  in  the  width  of  the  directivity  pattern  must  be  ex- 
pected.  It  is  also  noted  that  the  low  frequency  component  is 
sensitive  primarily  to  the  spacing  between  the  two  elements  of 
each  pair.   The  spacing  between  different  product  pairs  in  the 
array  influences  the  sidelobe  level  through  the  (sin  nBx)/(nBx) 
coefficients.   For  example,  consider  the  two-product-pair  cor- 
relation array 


'         *  3 

V«—  dlt  —*7^ —  4*s  — *T*"~  **** 

L^j> 1  I ^x} 

Ltfoj *TxU luTftsa 

—      Jl     


—  JTsin  i 

Vout  ^ti    XI2B 

[sin  Xl 
LX123 


T 


I — *{*> — I 


S3 


sin   X      B        sin  X    ,B 
+  JJt_ 


X34B 


X14B 


Sln  Xll!jcos    (X12+X34)2f0   + 


X--j  nB 


nB      sin  X_.B        sin  X.  _B      sin  Xoy  b1  .  Kt%r 

+  13     .  24    [cos    (X12-X34)2f 


X34B 


*HB 


An  -y  3 


I 


with  X^j  =  (trd^j/c)  sin  ^.   Since  spacing  between  the  pairs 
has  only  a  secondary  effect  on  the  directivity  pattern,  this 
can  be  eliminated,  giving  two  product  pairs  from  three  elements 
The  directivity  pattern 
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Ibofosh 


LoPgS5 


5 

j 


then   given    by    (with  relative    spacing   D   and    3D): 

^Jfosin   XB   sin   3XB|rr,o  j(tf  v\   + 

^^Tt^L  "TS ITT- JCOS  2C  0 ' 

[sin  XB.sin  gXB  +  sin  4XB|cos  4(2fnX)j  . 
XB      3XB      4XB    J  J 

As  a  comparison  of  the  directivity  patterns  which  re- 
sult from  the  general  correlation  array  utilizing  all  fre- 
quencies and  from  the  low  frequency  correlation  array,  consider 
a  four-element  array  under  both  conditions: 
■       Z  3 


■        Z  3  * 

HSJ* — 1 H3*—J K3« — J 


I 45k 


2 


With  relative  spacings  d^  =  4D,  d23  =  D,  d-,  m   2D,  the  equations 


for  the  two  cases  are 
1.  General  case 


out 


cc[A  cos   9(2fQX)   ♦  4B  cos   7(2fQX)   +  C  cos   5(2fQX)   + 


2D  cos   3(2fnX)   +  2E  cos    (2fnX)]  , 


23. 
2.  Low  frequency  case 


Vout   oc  [F   cos    7(2f0X)    +   G   cos    5(2fQX)    +  H  cos   3(2fQX)    + 

I    cos    (2fQX)] 

In  each  case  the  coefficients  A,  B,  ....  are  functions  of  (sin 
n  XB)/(n  XB)  whose  values  are  approximately  unity  for  narrow 
band  situations  and  become  less  than  one  only  in  wide-band  cir- 
cuits.  Assuming  a  value  of  unity  for  these  coefficients,  the 
directivity  patterns  for  these  two  cases  are  shown  in  Fig.  3. 
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25. 

III.   EFFECT  OF  FINITE  AVERAGING  TIME 

The  antenna  directivity  patterns  in  correlation  arrays 
are  exact  only  in  the  limit  of  infinite  averaging  time.   In 
actual  correlators  this  integration  time  would  be  a  finite 
interval  chosen  as  a  compromise  between  the  desired  response 
of  the  system  and  the  maximum  time  delay  which  can  be  tolerated 
It  becomes  important,  then,  to  investigate  the  system  per- 
formance as  a  function  of  this  integration  time. 

An  expression  for  the  effect  of  an  averaging  device  can 

be  obtained  from  general  filter  considerations. 9   If  h(t)  is 

the  effective  weighting  function  of  a  linear  measuring  device, 

and  x(t)  is  the  function  to  be  measured,  a  measurement  M  (T) 

x 

made  at  time  t  =  T  after  x(t)  has  been  introduced  at  t  =  0  can 
be  expressed  by  the  convolution 


T 
M  (T)  =  \h  (u 


)x(T-u)  du  h(u)  =  0,  u  <  0 

0 

where  (0,T)  is  the  observation  interval.   M  (T)  will  vary  from 
observation  to  observation,  fluctuating  about  the  expected 


value  M  (T)  with  a  variance  0"  =  OYF  -  [M  (T)]  . 
x  *        x 

In  the  general  situation,  x(t)  is  at  least  wide  sense 
stationary,  and 

cT 

F£Tr7  =  \  h(u)  x(T-u)  du 
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T  T 


M^CT)   =  \  V  h(u)  x(T-u)xCT-v)  h(v)  du  dv 


0  0 
T  T 


=  \  \  h  (u)  yx(a-v)  h(v)  du  dv 


0  0 


where  *f>  is  the  correlation  coefficient  of  the  function  to  be 
measured.   For  an  ideal  integrator  h(u)  =  1/T,  0<  u<  T.  And 

mtft  ■  ^rry 


T 


0 

This  result  was  obtained  by  Davenport*-"  in  an  early 
study  of  correlation  detectors.   However,  his  analysis  was 
directed  to  a  description  of  signal-to-noise  ratio  in  these 
detectors,  while  the  present  investigation  is  concerned  with 
placing  confidence  intervals  around  the  expected  value  of  the 
output  voltage,  so  the  methods  of  utilizing  the  information 
in  these  equations  will  differ. 

The  information  required  for  this  statistical  descrip- 
tion of  the  output  voltage  will  be  the  expected  value  and  the 
autocorrelation  coefficient  of  the  input  voltage  to  the  time 
averager.   The  calculations  will  be  discussed  in  some  detail 
for  the  two-element  array,  and  then  the  effect  of  more  elements 
in  the  correlation  array  will  be  indicated. 
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If  we  consider  the  usual  two-element  array,  an  arbi- 
trarily located  signal  source  induces  a  band- limited  random 
signal  V,(t)  on  element  #1  and  V  (t)  »  V.(t  +T)  on  element 
#2.   Then,  after  multiplication 

Y(t)  -  V1(t)V2(t)  =  V1(t)V1(t  +T) 

7TE7  -  R  (X) 

Vv(s)  =  Y(tWt+s)  -  V1(t)\/I(t+T>VI(t+s)V1(t4*+s5 
=  \(X)2   +  ^(s)2  ♦  RyCK*  s)Rv(1f-  s) 

The  last  term  of  this  expression  for  the  autocorrelation 
coefficient  cf  the  product  voltage  becomes  rather  awkward  to 
handle  in  the  equations  for  the  integrator  circuit.   And  as 
more  elements  are  added  to  the  array,  the  corresponding  terms 
in  their  product  equations  increase  the  complexity  of  the 
calculations.   Since  the  interest  here  is  in  the  change  in  in- 
tegration time  required  as  the  number  of  elements  in  the 
correlation  array  increases,  an  approximation  for  the  corela- 
tion  coefficient  will  be  used  which  will  permit  straightfor- 
ward calculation  of  the  integration  characteristics. 

If  the  signal  source  is  located  in  the  principal  lobe 
(t  ■  0),  the  equations  for  the  product  voltage  become 


Ytt7  -  Rv(0) 

Y„(s)  -  R(0)2  ♦  2R(s)2 

r 


=  R(0)2(1  ♦  2/Ks)2) 
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where/)  (s)  is  the  normalized  autocorrelation  coefficient; 
/O(s)=R(s)/R(0) .   Then  for  an  arbitrarily  located  source,  the 
approximate  product  voltage  equations  are 

yTT7  =  yr) 

¥>v(s)»Rv(<s)2(l  +  2/o(s)2) 

This  approximation  is  good  for  positions  near  the  principal 
lobe  axis  and  fails  for  positions  that  make  Ry(T)  quite  small 

(Rv(r)-*o). 

For   the  usual   rectangular   band-limited  voltage, 
>Ks)   =    (sinirBs/trBs)   cos   2irfQS,    and   the   equations   for  the 
ideal   integrator  are 


M^Ct)   =  RV(T) 


M  (T)      •  R   (X)2  f  \    (1   -    S)A+2fiin2n|s   cos22Tf0s)   ds 
T 

1  c 

my(t)    =  Rv(t)2  |  \  (i  -  |)[i+[^aafiflj2  * 


/sin*Bs\ 


2 


cos  2ir(2f0)s]ds 


The  solution  of  this  integral  can  be  expressed  in  two  parts, 
the  first  depending  on  the  low  frequency  components  of  the 
product  voltage,  the  second  determined  by  the  doutie  frequency 
component . 
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MY 


(T)2   =  RvCT)2[l   +  mj   +  m22] 


Or,  in  terms  of  the  variance  at  the  output  of  the  integrating 
circuit, 

<T2   =   R^T)2    [mj   ♦  m22] 


where 

m   2 
If 


=        .    1    -    cos    2»BT     +  JLSi(2lr3T) l-*Cin<2irBT>l 

L  (»BTr  WBT  (-BT)2  ] 

^~  ,for   BT  »    1. 
m  2   =  fS1(^5T>        sl(2kwBT?        <k+l)Si<;2k+2>rBT     + 

2f     I     Vbt  tTET  SFbt 

(fc-OSl»k-2>BT  +  sinW  cos2k.BT  .  Cin(2»BT)  _ 
2*BT  ^.BT)^  2(irBT)2 

Cin(2k»BT)  +  Cln(2k+s)TiBT  +  Cin(2k-s)r3T*] 
2(iTBT)2       4(WBT)2        4(fTBT)2    J 

where  k  ■  (4f  )/B,  Si(x)  is  the  sine  integral,  and  Cin(x)= 
Vl-cosv/v)  dv. 

The  standard  deviation  of  the  voltage  at  the  output 
of  the  integration  circuit  as  a  function  of  the  integration 

time,  T,  is  shown  in  Fig.  4.   In  the  limit  as  T— »0,  the 

2        .2 
variance  is  2  R  (t)   =  2M  (O  ,  as  would  be  expected  for 

the  Chi-square  distribution.   For  voltage  samples  observed 

without  integration,  a  95  per  cent  confidence  interval  on 

the  output  voltage  extends  from  zero  to  3.87  M  (T) .   As 
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integration  time  increases,  the  confidence  interval  becomes 
smaller,  and  observed  voltage  samples  would  cluster  more  and 
more  closely  about  the  expected  value,  M  (T).   Figure  5  shows 
this  narrowing  of  the  95  per  cent  confidence  interval  as 
integration  time  increases. 

It  should  be  noted  that  the  double  frequency  component 
of  variance  decays  rapidly  compared  to  the  low  frequency  com- 
ponent.  A  correlator  whose  integration  time  is  long  compared 
to  1/f  may  still  have  an  output  voltage  with  an  appreciable 
variation  because  of  this  low  frequency  variance.   If  the  input 
circuits  of  the  integrator  rejected  this  double  frequency 
component  of  the  product  voltage,  the  output  voltage  would  fluc> 
tuate  with  only  the  low  frequency  component  of  variance. 
Regardless  of  this  circuitry,  however,  the  integration  time 
that  will  be  necessary  to  reduce  the  fluctuation  will  be  very 

long  compared  to  1/f  ,  and  the  double  frequency  component  of 

0 

variance  can  generally  be  completely  neglected. 

Incidentally,  this  statistical  description  of  the  fluc- 
tuation of  the  output  voltage  of  a  two- element  correlation 
array  is  the  same  as  the  statistical  description  of  the  signal 
power  output  of  a  linear  additive  array.   The  only  difference 

would  occur  in  the  magnitude  of  the  expected  value,  M  (T) . 

Y 
Knowing  the  variance  of  the  linear  array  power  fluctuation  as 

a  function  of  time,  we  can  compare  averaging  time  requirements 

of  more  complex  correlation  arrays  with  that  of  the  familiar 

additive  array. 
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The  autocorrelation  coefficient  of  the  product  voltage, 
Y(t),  for  general  correlation  arrays  will  be  calculated  in  the 
same  manner  as  for  the  two- element  array.   The  expected  value 
of  the  product  voltage  can  be  determined  from  the  charcteristic 
function  as  in  the  preceding  section.   The  autocorrelation 
coefficient  of  the  product  voltage  for  a  source  located  on  the 
main  lobe  axis  can  also  be  determined  by  the  charcteristic 
function  mdhod.   This  will  be  modified  by  inserting  the  ex- 
pected value  of  the  arbitrarily  located  source,  and  the  result- 
ing form  will  be  used  to  approximate  the  actual  autocorrela- 
tion coefficient. 

The  expected  value  of  product  voltage  is  given  in  the 
following  table  for  a  few  values  of  the  number  of  elements,  n. 

2  R(t  ) 

12 

4     R(T10)R(r/)+R(r  )R(r  )+r(t  )R(r  ) 

12    34     13    24     14    23 
6 

6      ^T    R(T  )R«  )R(*  ),  ra*n#oj«p*q 
r  lm    no    pq 

m,n,o,p,q=r2 
The  characteristic  function  for  the  product  voltage  i6 
F<f  i.f2»  s)  *  E[>exP   CJ  /tVCt:)  +  j/2V(t+.)]J 

Then  Vn(s)  is  the  coefficient  of  (jf  ^"/n  -  (j/2>n/n.  in 
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the  power  series  expansion  of  F(£ i  >  f7  >  s). 

n  y  (s) 

2  VTtT2  [1  +  2  *(s)2] 

4  YTFT2  [9  +  72  *(s)2  +  24*(s)4] 

6  YTtT2  (225  +  4050^0  (s)2  +  5400  P(s)4  ♦  720*(s)6] 

Calculation  of  the  variance,  for  the  usual  rectangular  band- 
limited  voltage,  gives 

2   7  l 

for  two  elements  (and  for  signal    f     -   Y(t)   -i- 

power  in  linear  array)  Y  Bt 

for  four  elements  j-  =  Y(t)Z   •%— 


Y 


BT 


for  six  elements  *-  =  Y(t)       ' 

Y  BT 


in  all  cases,  BT  »  1 

This  last  table  shows  the  general  relationship  between 
variance  and  number  of  elements  in  a  correlation  array.   For 
any  established  level  of  fluctuation  in  the  output  voltage,  a 
six-element  correlation  array  will  require  an  integration  time 
15.5  times  as  long  as  is  required  for  a  linear  array  (or  a  two- 
element  correlation  array).   In  the  previous  section  the  direc- 
tivity pattern  of  a  correlation  array  was  found  to  be  similar 
to  that  of  a  much  longer  linear  array.   Thus  the  correlation 
process  in  antennas  is  seen  to  be  essentially  one  of  trading 
space  for  time. 
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IV.   RESOLUTION  CHARACTERISTICS  OF  CORRELATION  ARRAYS 

In  describing  the  capability  of  a  correlation  array  to 
resolve  two  similar  signal  sources  separated  by  an  angular  dis- 
placement 0,  the  directivity  patterns  developed  in  Section  III 
cannot  be  applied  directly.   Since  the  correlation  arrays  de- 
pend on  multiplicative  operations,  the  presence  of  two  sources 
causes  cross-product  terms  to  appear.   These  terms  then  neces- 
sitate additional  calculations  to  determine  resolution  capa- 
bility.  These  calculations  will  be  examined  for  correlation 
arrays  in  the  presence  of  single  frequency  sources  and  randomly 
varying  sources.   The  correlation  arrays  will  be  found  to  have 
resolution  capabilities  equivalent  to  those  of  longer  linear 
arrays,  if  the  sources  to  be  resolved  are  independent;  however, 
for  coherent  sources  resolution  may  be  possible  only  for  cer- 
tain specific  separations  of  the  sources  and  not  as  a  general 
rule. 

First,  let  us  examine  a  four-element  correlation  array 
as  developed  in  Section  III.   Assume  there  are  two  raonochromatic 
sources,  source  A  (V  ■  A  cos  u*xt^  anc*  source  B  (V*  a  cos««»  t), 

separated  by  an  angular  displacement  0  .   The  general  expression 
for  the  voltage  response  of  the  array  as  a  function  of  position 
of  the  sources  is  quite  lengthy,  and  is  developed  in  Appendix 
II.   It  is  found  to  depend  on  the  coherence  of  the  sources. 
The  two  cases  that  occur  can  be  summarized: 
1 .  For  two  coherent  sources : 
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4  / 

VQut   oc    A    (cos    X   +   cos    3X   +    cos    5X)    +    B4(cos    Y  +   cos    3Y  +   cos    5Y) 

+A3B    £  [cos(nX-raY)]    +    A2B2f [cos(nX-mY)]    +  AB   f  [cos(nX-mY)  ] 

2.   For  two  sources  of  slightly  different  frequencies: 

VoutocA4(cos   X  +   cos    3X  +   cos    5X)    +    B4(cos    Y  +   cos    3Y  +   cos    5Y) 

+A3B   f  [cos(nX-mY)-    cosSt]*  AB3f [cos(nX-mY) • cos  St] 

2    2  7    7 

+A   B   f [cos(nX-mY)-cos  it]    +  A.  B   f  [cos    (nX-mY)] 

where  X  =  (u»D/c)  sin  ^A  and  Y  =  (a>D/c)  sin  ^   are  the  phase 

delay  factors  for  the  two  sources  with  element  spacing  D,  and 
f  [cos (nX-mY) ]  represents  a  number  of  cosine  terms  containing 
these  phase  delay  factors,  while  f[  cos  (nX-mY)*  cos  St]  repre- 
sents a  number  of  these  cosine  terms  modified  by  the  beat  fre- 
quency component,  St,  between  source  A  and  source  B. 

In  both  of  these  equations,  the  first  two  terras  corres- 
pond to  the  voltage  response  pattern  of  a  linear  antenna  array 
in  the  presence  of  two  sources.   The  remaining  terms  arise  from 
the  nonlinearity  of  the  correlation  array.   In  the  ca9e  of  two 
sources  with  identical  frequencies,  the  cross-product  terms 
are  constant  with  time  and  averaging  the  output  roltage  will 
not  alter  the  result.   When  the  sources  have  different  fre- 
quencies, however,  time  averaging  can  be  employed  to  reduce  the 
beat  frequency  (cos  St)  part  of  the  cross-product  terms;  but 
even  in  this  case  the  resulting  voltage  response  pattern  will 
not  be  precisely  equivalent  to  that  of  the  linear  array,  since 
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some  of  the  cross-product  terms  will  remain  constant  with  time 
and  will  contribute  to  the  dc  component  of  the  output  voltage. 

The  voltage  response  pattern  of  this  four-element  cor- 
relation array  as  the  array  is  rotated  past  a  pair  of  coherent 
signal  sources  separated  by  an  angular  displacement  of  20°  is 
shown  in  Fig.  6.   In  this  case  the  cross-products  result  in  an 
apparent  indication  of  a  single  source  located  midway  between 
the  actual  pair  of  sources.   This  four-element  correlation  array 
will  resolve  sources  with  an  angular  separation  of  14°  and  with 
an  angular  separation  of  28°,  but  a  pair  of  sources  whose  angu- 
lar separation  is  in  the  range  0°-13°  or  in  the  range  15°-26° 
will  not  be  resolved  by  the  array. 

It  should  be  noted  that  if  the  correlator  contains  only 
one  stage  of  multiplication,  then  the  cross-product  terms  oc- 
curring from  two  sources  of  different  frequencies  appear  only 
as  beat  frequency  components.   These  terms  describe  a  low  fre- 
quency ac  signal  appearing  with  the  desired  dc  measurement 
voltage.   In  this  case  the  cross-product  terms  can  be  minimized 
by  time  averaging  the  output  voltage.   However,  if  more  than 
one  stage  of  multiplication  occurs  between  the  antenna  element 
and  the  output  of  the  array  (as  is  the  case  in  the  example  of 
the  four-element  correlation  array),  the  cross-product  terms 
will  occur  as  low  frequency  ac  terms  and  also  as  dc  terms.   In 
this  case  the  equivalence  between  the  correlation  arrays  and 
the  linear  arrays  is  not  immediately  apparent  but  must  be 
determined  by  calculation  of  resolution  characteristics  for  the 
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particular  type  of  signal  encountered. 

Assuming  a  basic  spacing  D  =  X/2,  the  four-element  cor- 
relation array  of  length  4D  will  resolve  two  sources  of  slightly 
different  frequencies  at  approximately  19.5°.   This  is  equi- 
valent to  the  resolution  of  a  uniform  linear  array  vith  an 
aperture  of  9D,  about  twice  as  long. 

A  slightly  more  complex  example  of  a  correlation  array 
is  one  discussed  by  Drane : 


I 


J 


The  element  on  the  left  represents  a  uniform  linear  array  with 
aperture  2a,  while  that  on  the  right  is  a  simple  interferometer 
with  aperture  2b. 

The  directivity  pattern  of  the  uniform  linear  array  is 
proportional  to  (sin  X)/X,  where  X  =*>(a/c)sin  tfl.  The  direc- 
tivity pattern  of  the  interferometer  is  proportional  to  cos  Y, 
where  Y  =*^(b/c)sin  4.  If  a=b  and  there  is  no  spacing  between 
the  right  end  of  the  linear  array  and  the  left  element  of  the 
interferometer,  the  directivity  pattern  of  the  correlation 
array  is 

«    «.  sin  4X 
V     cc    . 
out     4X 


which  is  the  same  as  the  directivity  pattern  of  a  uniform  linear 
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array  of  twice  the  length  of  this  correlation  array. 

However,  if  two  coherent  sources  are  present  the  output 
becomes 

sin  4X    sin  4X  sin  x  x  cos  x  sin  x 

Vout  oc. *  ♦ 2  +  cos  (X.  +  X  )      A     Bt    _A b 

4X.       4X  ah     4X  4X 

A3  A  B 


The  third  term  again  occurs  because  of  the  nonlinearity  of  the 
correlator.   Calculation  of  the  resolution  capability  of  this 
array  shows  that  it  is  equivalent  to  a  uniform  linear  array 
1-1/2  times  as  long  as  the  correlation  array. 

These  have  been  only  two  examples  of  the  effect  of  cross- 
product  terms  in  a  cotrelation  array  in  the  presence  of  mono- 
chromatic sources.   In  each  case  that  might  t>e  considered,  the 
effect  of  these  terms  would  vary,  depending  on  the  types  of 
signal  voltages  emitted  by  the  sources  and  by  the  number  of 
successive  multiplicative  processes  in  the  correlator. 

When  two  statistically  similar  sources  are  inducing 
randomly  varying  voltages  on  the  antenna  elements,  the  reso- 
lution capabilities  of  the  correlation  array  must  be  expressed 
in  terms  of  integration  times  and  confidence  intervals.   Con- 
sider again  the  four-element  correlation  array  and  two  inde- 
pendent sources  (source  A  and  source  B).   Then  the  expected 
value  of  the  output  voltage  can  be  calculated  either  directly 
or  by  the  characteristic  function  method  and  can  be  expressed 
in  terms  of  the  correlation  coefficients  of  the  individual 
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voltages.   The  voltage  induced  on  the  i-th  element  is 

V.  (t)  =  V.(t  +  T.)  +  VR(t  +  X   ) 
and 

vTt7  -  VL(t)V  (t)V3(t)V  (c)  = 

+    [RB(tl2>RB(T34)+  1B<'«i3>Vr24>4  RB<'«  u)RB<'e23)] 

♦  [RA(x12)RB(x34)+  »A( *„)*,<  V*  W*8(X235 

+     V*34)RB(T12)+  WW  RA(T23)V^4) 

The  first  and  second  lines  in  this  expression  give  the  expected 
voltages  due  to  each  source  individually.  The  third  and  fourth 
lines  contain  the  cross-product  components  which  occur  as  a  re- 
sult of  the  two  stages  of  multiplication  of  the  element  voltages 

The  variance  of  this  output  voltage  could  be  calculated 
directly;   however,  the  large  number  of  terms  in  the  final  ex- 
pression would  make  this  quite  laborious.   Therefore,  this 
quantity  will  be  approximated  as  previously  by  considering  the 
variance  of  the  voltage  produced  by  a  single  source  on  the  main 
lobe  axis  and  then  introducing  the  new  value  of  Y(t )  given  by 
the  equation  above. 

With  an  arbitrary  angular  displacement  between  the  two 
sources,  and  a  given  level  of  confidence  desired  in  the  result, 
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it  is  possible  to  calculate  the  maximum  variance  which  will 
permit  a  significant  difference  in  the  amplitude  of  the  output 
voltage  when  the  array  is  directed  at  a  source  and  when  it  is 
directed  between  the  sources.  This  value  of  variance  can  then 
be  used  to  determine  the  minimum  permissible  integration  time. 
So  the  array  resolution  is  a  function  of  two  quantities:  level 
of  confidence  and  integration  time. 

The  resolution  capabilities  of  the  simple  four-element 
correlation  array  are  listed  in  the  following  table.   They 
are  compared  with  those  of  a  linear  array. 

Random  Signals 
Resolution  at  the  95  Per  Cent  Confidence  Level 


Four- Element  Correlation  Array 

Length  4D 

Resolution  Integration  Time 

19°  BT  =   18 

18°  BT  =   97 

17.5°  BT  =   890 

17°  BT  =   od 


Aperture  of  Equivalent 
Uniform  Linear  Array 


9.6  D 
9.4  D 
9.2  D 
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V.  EFFECT  OF  DISTRIBUTED  NOISE  ON  THE  NONLINEAR  ANTENNA  SYSTEM 


The  following  assumptions  will  be  made: 

1.  uniform  distribution  of  independent  random  noise 
sources; 

2.  noise  sources  in  an  element  d-H-  induce  a  noise  power 
#-«*•  watts /cps  in  an  isotropic  antenna  element; 

3.  antenna  element  and  associated  circuitry  have  a  rec- 


1 


tangular  passband  of  bandwidth  B,  center  frequency  f  ,  with 
fQ»  B. 

Then  the  autocorrelation  function  of  the  noise  voltage 
induced  on  an  antenna  element  is 


Vu   (s)  =  xL(t)x1(t+s)  =  4i»  Bn  "7b 

=  *n 


sinrr  Bs 


coso)0  s 


(0)  sin  n  gs  cos  cO.  s 
ttBs 


where  R^O)  represents  the  average  noise  power. 

With  two  identical  isotropic  elements,  separated  by  a 
distance  d,  the  voltage  induced  on  one  element  by  the  noise 
sources  in  an  angular  area  d-fl-  differs  from  that  induced  on 
the  second  antenna  element  only  by  a  delay  time  X 


X  =  —  cos  e 
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The  same  average  power,  R^CO),  is  induced  on  each  element.   The 
autocorrelation  coefficient  of  the  noise  voltage  on  each  element 
is  the  same. 

The  cross-correlation  coefficient  becomes: 

Vl2(s)   =  #1(t)x2(t+s)  =  x^tjx^t+t  +s) 

-  ffit'^'ftff  cosul.(r+s)  sin  0  dS  drf 

=  2wB»(^ln,rB(c  cosO+s) 

ll*±   cos  0+s)   C°S  ^(c  COsG+8)  »in  •  d  0 

This  can  be  integrated  to  yield: 

Y\2<s>  =  K(0)^-(Si[(k+1)Ei1-+irBs)]  -Si[(k+l)(a|L+lBs3 

r(k-i)(J^  +rBs)l    +  Si[(k-l)(zlBd  +7r  Bs)j[ 

where  k  =  <***     and  Si(x)  is  the  sine  integral. 
1TB 

By  the  symmetry  of  the  integrand  and  of  the  range  of 

interpretation  it  is  evident  that 

V21(S)  -  X2(t)x1(t+8)  -  XL(t+*)x1(t  +  8>  *^12Cs). 


Under  the  usual  conditions  of  f  »B,  and  for  other  than 

o 

very  close  element  spacings,  this  expression  for  the  cross- 
correlation  coefficient  can  be  considerably  simplified.   Under 
these  conditions,  approximate  Si(x)  by  TJ  ■  ffl8  x,  and,  after 


-Si 
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some  manipulation,    find 


<f12(s)   =  RnCO)    sin  kJ^  ^^        sinirBs 


]*C|^) 


TTBd 


FBI 


COS  44,  s 


where 


*  Rn<t)     linirBs     cos  s 

ttBs 


R(t)  =R(0)  sltVsVL) 


n 


n 


oyd 


L°feM) 

T-Bd 
c 


Since   in  most   cases   k»  1,   while  d,    at   the  most   is   a   few  wave- 
lengths,   the   last   factor,    sinC*  Bd)/    U^&)  s   win   have  very 

c 

little  influence  on  the  expression;  and  the  cross-correlation 
coefficient  can  be  considered  to  be  a  function  of  sin(^*—  )/(-^r) . 

It  should  be  remembered  that  the  approximation  used  to  simplify 
the  sine  integrals  in  the  original  equation  is  not  valid  for 
small  values  of  element  spacing,  d;  and,  in  the  event  of  close 
element  spacings,  calculation  of  the  cross-correlation  coefficient 
would  require  use  of  the  exact  equation. 

Calculation  of  noise  voltage  and  noise  power  present  at 
the  output  of  a  simple  nonlinear  array  will  demonstrate  the 
considerations  which  must  enter  into  a  specification  of  antenna 
element  spacing.   Take  for  example  a  four-element  array: 


I 


3Y 


n»<t) 


•  x 


-E- 


*.«> 


Tver) 
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The  noise  sources  are  assumed  Co  be  uniformly  distributed  over 
the  surface  of  a  unit  sphere  surrounding  the  array.   Then  the 
average  value  of  output  voltage,  Y ( t ) ,  becomes: 

YTFT  =  VTtT  =  n1(t)n2(t)n3(t)n4(t). 

=  VT12)Rn(<*34)  ♦  ^n^l3)Rn(r24)  +  Rn(  ^  14)R-(^23  > 

Then,  as  shown  in  a  previous  section,  the  average  value  of  out- 
put power  is: 

T 


^n- 


Y(t?  =  i    \     (1  -  f )  V(t)v(t+s)   ds 


But,  with  four  elements, 

V(t)V(t+s)  =  ni(t)n2(t)n3(t)n4(t)n1(t+s)n2(t+s)n3(t+s)n  (t+s) 

And,  when  expressed  as  multiples  of  the  correlation  coefficients 
of  the  element  noise  voltages,  this  results  in  an  equation  with 
some  105  terms.   However,  if  spacings  are  selected  to  make 


H^jCs)  =  0,  (i  *   j),  then 


V(t)V(t+s)  -[Yn(s)]4  «  \(0)4(^p) 


4    4  fc 
cos  co»  s 


and  the  output  power  can  be  readily  evaluated.   Spacings  which 
give  this  result  are  multiples  of  a  half-wavelength  at  the  center 
frequency,  f  . 

Since  the  cross-correlation  coefficients  can  assume  both 
positive  and  negative  values,  it  is  conceivable  that  some  decrease 
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in  output  noise  power  could  be  realized  by  appropriate  changes 
in  element  spacings  after  a  study  of  the  complete  expression 
for  V(t)V(t+s)  under  the  existing  conditions  of  noise  source 
distribution.   However,  the  above  method  of  considering  all 
cross-correlation  coefficients  as  zero  gives  a  quick  evaluation 
of  output  power  which,  by  the  very  nature  of  the  sin  x/x  func- 
tion, is  subject  only  to  minor  modification  by  more  detailed 
considerations . 

This  development  of  the  output  noise  power  from  a  non- 
linear array  does  not  give  the  complete  description  of  the  oper- 
ation of  this  type  of  array  when  it  is  used  to  detect  a  signal 
buried  in  a  general  noisy  medium.   Because  of  the  multiplicative 
operations  performed  in  the  antenna  circuitry,  there  will  be 
an  interaction  between  signal  and  noise  which  can  appear  as  an 
undesired  component  of  the  output  voltage. 

This  interaction  will  be  shown  by  considering  a  simple 
two-element  correlation  array.   Assume  a  single  signal  source 
induces  a  voltage  which  is  in  phase  on  each  element,  and  a 
uniform  spherical  distribution  of  noise  sources  induces  unde- 


sired noise  voltages  on  the  elements. 

I 


2 
V 


V,(c) 


1  AV6.  | 


VtCT) 


Then 
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V]_(t)  =  S(t)  +  nL(t) 


V2(t)  =  S(t)  +  n2(t)  =  S(t)  +  n^t+T   ). 


YTtT  =  S2(t)  +  S(t)  nL(t)+  n  (t)   +  ni(t)n2(t) 


=  S2(t)  +  Rn(TL2)  =  S2(t)  +4,12(o) 

. 
If  the  element  spacing  is  chosen  (multiple  of  a  half  wavelength) 

to  make  4^l9(0)  =  0,  then  the  voltage  in  the  output  circuit 

due  to  the  signal  is 

Y(0  =  S2(t). 

Noise  power,  N,  in  the  output  circuit,  can  be  determined  from 

N  =  ±£  (1  -O  [Y(t)Y(t+s)  =  YTtT2]ds. 

TJo      (1   -  1)\   S2(t)S2(t+s)    +   S(t)S(t+s)  Tn^tJi^lt+sJ    + 

n2(t)n2(t+s)  I  +  nI(t)n2(t)nI(t+s)n2(t+s)    -S2(t)S2(t)V   ds 

retaining  only  terms  whose  expectations  are  nonzero. 

If  the  signal  voltage  induced  on  a  single  element  is 
a  simple  sinusoid,  S  =  A  cos  oj   t,  then  the  signal  voltage  in 
the  output  circuit  is 

YTtT  =  hk1. 

The  signal  power  in  the  output  circuit  is 

p  =  u\ 
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and   the   noise   power    is 


Vo)2J^iMcosU)|jds 

Since  f0»B,  the  ac  component  of  the  noise  will  be  essentially 

zero  in  averaging  times  such  that  BT  >  1 .   So  only  the  low 

frequency  component  of  the  noise  power  in  the  above  equation 

need  be  considered. 

T 
A       (1    -   |))AXC»  Sinj^s   +Rn(0)2      sln^Bs} 
T  t)  J ^ irBs  — j (irBs^j  ds 


N-l 


For  large  values  of  BT,  this  can  be  approximated  by 

N*  A2Rn(0)   1   +  ICCO)2  1_ 

- BT    -^ BT 


So  the  noise  power  appears  not  only  as  a  result  of  the  noise 
sources  alone  (last  term)  but  also  as  a  result  of  the  inter- 
action or  signal  and  noise  voltages  (middle  term). 

This  total  noise  power  resulting  from  the  interaction 
of  signal  and  noise  voltages  may  well  be  greater  than  that 
caused  by  the  noise  sources  alone.   For  this  two-element  non- 
linear array  and  under  conditions  of  unity  signal-to-noise  ratio 
on  an  individual  element,  the  total  noise  power  resulting  from 
the  cross  multiplication  is  twice  the  noise  power  caused  by 
the  noise  sources  alone. 

Nonlinear  arrays  of  this  same  type  (combining  network 
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composed  of  product  circuits  only),  but  with  greater  numbers 
of  elements,  will  have  additional  noise  power  terms  resulting 
from  this  cross  multiplication  of  signal  and  noise  voltages. 
This  is  summarized  in  the  following  table.   In  this  table, 
the  level  of  multiplication  refers  to  the  number  of  times 
the  voltage  induced  on  an  element  undergoes  multiplication  be- 
fore reaching  the  output  circuit. 
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VI.   NONLINEAR  ANTENNAS  AND  CHANNEL  CAPACITY 

A  direct  comparison  of  the  effect  on  channel  capacity 
of  the  linear  and  nonlinear  antennas  is  complicated  by  the 
change  in  dimensionality  of  the  message  space  effected  by  the 
multiplicative  and  averaging  processes.   A  comparison  is 
possible,  however,  between  the  nonlinear  antenna  system  and 
a  linear  antenna  with  an  associated  square  law  circuit.   The 
signal-to-noise  ratios  in  the  outputs  of  these  two  systems 
can  show  the  effectiveness  of  the  nonlinear  antenna  system 
as  related  to  the  familiar  linear  antenna  system. 

Consider  an  antenna  array  of  p  elements,  spaced  so  the 
random  noise  voltages  induced  on  any  two  elements  are  un cor- 
related.  Then,  if  s  is  the  voltage  induced  on  each  element 
by  a  signal  source  in  the  main  lobe,  and  n.  is  the  random 
noise  voltage  induced  on  the  i-th  element,  the  voltage  from 
the  linear  array  is 

P 
V,  *  ps  +  2l2  n 
L       i  =  1  i 

After  squaring  and  averaging 

Y  «  V2  =  p2"?  ♦  pRn(0). 
L    L  n 

But  if  these  p  elements  are  split  into  two  groups  of  p.  and  p 

elements  respectively  (p^+p2  -  p)»  *n<*  the  voltages  of  these 

Cwo  sub-arrays  are  multiplied  and  averaged,  the  output  of  this 

nonlinear  array  becomes: 


■  (p,  s  +  ^~7  n-  )(p 
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— _ 

_2 


NL    "1     i=l 


. )(p0s  +  SET  n  _ )  =  p  p  s 
l   2    j  =  i   j      J-  I 


The  noise  power  in  the  output  ot  each  system  will  depend 
on  the  type  of  signal  source.   Two  such  sources  will  be  examined 
(.1)  a  random  signal  source  with  a  g^u     .  .  .str  .but ion,  ana 
{L)    a   simple  sinusoid,  •  =  A  cos   t.   Then  in  terms  ot  the 

signal,  to  noise  power  ratio,  £, 

N 


1.   s  gaussian. 

Linear  Array  u>)    =  gT 
N  out 


(i}  in  \(5)2  +  2  sx TT-  \ 

L&  i"  p^in   p2  J 


Nonlinear  Array  (£)   =  BT(§)2  \  2PlP2 


out     iN   in 


2PlP2(|)^n  +  p(|)  +  1 


when  p^  =  P2  =  P 


(!>out  ■  »<i>2i»  [  (|>L  ♦  fcf>  ♦  J 


2.   s  =  A  cos  w  t. 


Linear  Array  (S)   .  =  BT(S)2    T — 7-7-5 

N  out       N  m  I  — +   +  ^(-V 

L  P2   pvN'in 

U  ♦*(!>,.    J 

LD    p  N  in  J 


Nonlinear  Array  (§)out  =  B^ 

when  Pi  =  Po  =  § 

(S)    =  BT(|)2 
N  out      N  Ln 


where  (£)    refers  to  the  signal-to-noise  power  ratio  in- 
N  in 

q 

duced  on  a   single   element  of   the   array,    and    (-)        refers   to 
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the  signal-to-noise  power  ratio  at  the  output  of  the  antenna 

array  combining  circuit. 

When  the  signal-to -noise  ratio  and  the  number  of  elements 

are  such  that  their  product  is  large  I  p(^)   >>   ll  ,  the  out- 

L   N  in      J 

put  signal-to-noise  ratio  is  approximately  the  same  for  both 
arrays : 


(§)    ^   3T(|)   (2\ 
N  out       N  in  V  2  / 


In  this  case,  the  channel  capacity  would  be  expected  to 
be  the  same  for  both  arrays  if  the  same  number  of  elements, 
p,  were  used  in  each. 

When  the  signal-to-noise  ratio  and  the  number  of  elements 
are  such  that  their  product  is  small   |p(^)^„<<  l|  ,  the  out- 
put signal-to-noise  ratios  become: 


[■4>in«    l] 


Linear   Array      <f >Qut  »     BT<f)in(p   ) 

9      9 

Nonlinear  Array  (f)^  ~     BT(§)ln(|) 

Now  the  channel  capacity  of  a  nonlinear  array  is  not  so  great 
as  that  of  a  linear  array  with  the  same  number  of  elements,  p. 

This  analysis  has  tended  to  ignore  one  very  important 
aspect  of  the  nonlinear  array.   The  average  output  voltage  of 
the  linear  array  (Y. )  will  contain  a  component  proportional 
to  the  signal  and  a  component  proportional  to  the  noise;  the 
average  output  voltage  of  the  nonlinear  array  (Y  )  has  no  com- 
ponent  which  is  independent  of  the  signal.   This  emphasized 
the  real  difference  between  the  linear  and  nonlinear  arrays; 
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the  nonlinear  array  will  generally  be  capable  of  a  smaller 
channel  capacity  than  the  linear  array,  but  in  the  final  out- 
put voltage  the  effect  of  the  random  noise  sources  can  be 
made  arbitrarily  small  in  the  nonlinear  array,  while  this  is 
not  possible  in  the  linear  array. 
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the  signal-to-noise  power  ratio  at  the  output  of  the  antenna 

array  combining  circuit. 

When  the  signal-to  -noise  ratio  and  the  number  of  elements 

are  such  that  their  product  is  large  I  p(~).  »   l|  ,  the  out- 

L   N  in     -J 

put  signal-to-noise  ratio  is  approximately  the  same  for  both 

arrays : 

(f)    a*   BT(|)  (&\. 
^  out       N  in  V  2  / 

In  this  case,  the  channel  capacity  would  be  expected  to 
be  the  same  for  both  arrays  if  the  same  number  of  elements, 
p,  were  used  in  each. 

When  the  signal-to-noise  ratio  and  the  number  of  elements 
are  such  that  their  product  is  small   iP^m)-  <<a    H  >  tne  out- 
put signal-to-noise  ratios  become: 

9     9  v 

Linear  Array   (£)    *  BT(|) .   (p  ) 

J    ■    N  out       N  in  \   / 

9      9 

Nonlinear  Array  (f)^ »  BT(|)in(§) 

Now  the  channel  capacity  of  a  nonlinear  array  is  not  so  great 

as  that  of  a  linear  array  with  the  same  number  of  elements,  p. 

This  analysis  has  tended  to  ignore  one  very  important 

aspect  of  the  nonlinear  array.   The  average  output  voltage  of 

the  linear  array  (Y. )  will  contain  a  component  proportional 

to  the  signal  and  a  component  proportional  to  the  noise;  the 

average  output  voltage  of  the  nonlinear  array  (Y  )  has  no  com- 

NL 

ponent  which  is  independent  of  the  signal.   This  emphasized 
the  real  difference  between  the  linear  and  nonlinear  arrays; 
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the  nonlinear  array  will  generally  be  capable  of  a  smaller 
channel  capacity  than  the  linear  nrray,  but  in  the  final  out- 
put voltage  the  effect  of  the  random  noise  sources  can  be 
made  arbitrarily  small  in  the  nonlinear  array,  while  this  is 
not  possible  in  the  linear  array. 
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VII.   DESIGN  OF  NONLINEAR  ARRAYS 

A.   Basic  Considerations 
The  preceecling  analyses  form  the  basis  for  certain  con- 
clusions regarding  the  employment  of  antenna  arrays  which 
utilize  nonlinear  processes  in  the  combining  network: 

1.  Random  noise  voltages  induced  on  the  antenna  ele- 
ments can  be  reduced  to  any  desired  level  by  proper  element 
spacing  and  by  averaging  for  a  sufficiently  long  interval. 
Noise  power  is  approximately  inversely  proportional  to  the 
averaging  interval,  BT. 

2.  Cross  multiplication  processes  in  the  antenna  circuitry 
will  introduce  additional  noise  terms  resulting  from  inter- 
action of  signal  and  noise  voltages,  of  different  noise  voltages, 
or  of  different  signal  voltages.   These  added  noise  terms 

attain  their  minimum  values  when  the  level  of  multiplication 
is  limited  to  one. 

3.  When  the  level  of  multiplication  is  greater  than 
one,  the  interaction  of  two  or  more  independent  signal  sources 
will  introduce  components  in  the  output  voltage  which  can  not 
be  reduced  by  averaging. 

4.  The  examples  of  nonlinear  antennas  considered  have 
demonstrated  directional  patterns  with  main  lobes  significantly 
narrower  than  those  obtainable  from  linear  arrays  of  the  same 
size. 
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5.   The  nonlinear  array  is  not  efficient  in  a  general 
communication  system,  but  offers  particular  advantages  to  a 
system  in  which  narrow  beamwidth  is  paramount. 

These  general  conclusions  will  be  the  basis  for  the 
deisgn  of  nonlinear  antenna  arrays.   The  difficulties  which 
arise  when  the  level  of  multiplication  exceeds  one  outweigh 
any  apparent  increase  in  directivity  which  might  be  obtained. 
Therefore,  the  combining  network  will  be  a  combination  of 
addition  and  multiplication  circuits,  with  the  level  of 
multiplication  limited  to  one. 

B.   Products  of  Linear  Arrays. 

If  we  consider  a  uniform  array  of  n  equally  spaced  iso- 
tropic elements,  the  simplest  type  of  nonlinear  antenna  which 
can  be  constructed  (with  level  of  multiplication  equal  to  one) 
is  that  in  which  the  elements  are  formed  into  two  linear  arrays 
whose  output  voltages  are  then  multiplied  and  averaged.   There 
are,  of  course,  a  number  of  different  combinations  of  the 
elements  which  can  be  realized  in  accomplishing  this. 


V 


n,  element* 


58. 


For  two  groups,  containing  n   and  n  elements  respectively, 
so  that  n  =  n, +n  ,  and  with  spacing  d,  between  elements  and  d_ 
between  group  centers,  the  directional  patterns  can  be  written 


P   (»f/)  =  sin  ni(4V2) 
nl       ni  stn(f/i) 

Pn  (vp)  =  sin  n2(V^/2) 
2       n2  sin  Op/2) 

and    Pn(9)   =   sin  n^/2)    sin  n2(V/2)   .      , 

nL  sin(v/J)   n2  sin(y/2)   cos(£2y) 

The  two  parameters  controling  P  (<p)  are 

1.  the  total  number  of  elements,  and 

2.  the  number  of  elements  in  one  of  the  groups. 
Pn(M0  is  equal  to  unity  when  ty  =   Q .      Each  term  in  the 

expression  for  P  (4O  will  introduce  zeros  into  the  directional 

n 

pattern.   The  first  and  second  terms  will  haverzeros  equally 
spaced  along  the  tp  -axis  at  intervals  of  2TT/n^  and  llt/uy 
respectively.   There  are    (n^- l)*(n2-l)  =  n-2  I   such  zeros 
between  Y  =  0  and  4*  =  2TT  .   Assuming  there  is  no  overlapping 
of  the  two  groups,  d2/d^  =  n/2;  the  final  term  in  the 
directional  pattern,  cos(d2/d^)  ^  =  cos  nH^/2,  has  n  zeros  in 
the  interval  between  vp  =  0  and  9  =  2TT. 

It  is  evident  then  that  a  uniform  nonlinear  array  with 
n^+n^  =  n  elements  will  have  a  directional  pattern  with  2(n-l) 
zeros  instead  of  the  (n-1)  which  would  be  available  for  a 
simple  n-element  linear  array. 
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The  cos  n  *f/2  factor  in  the  directional  pattern  is  the 
important  innovation  of  this  type  of  nonlinear  antenna.   This 
factor  will  cause  the  pattern  to  pass  through  zero  for  V  =  T  /n, 
regardless  of  the  manner  of  splitting  the  n  elements  into  two 
groups.   If  these  n  elements  were  utilized  as  a  uniform  linear 
array,  the  first  zero  would  occur  at  f  =   2TT/n.   Thus  the 
beamwidth  of  the  principal  lobe  of  the  nonlinear  array  is 
approximately  \   that  of  a  linear  array  of  the  same  size. 

Since  the  location  of  n-2  zeros  will  depend  on  the  choice 
of  rij  and  n£  =  n  -  n.  ,  a  number  of  directional  patterns  are 
possible  for  any  fixed  total  number  of  elements,  n.   Figure  8 
compares  the  directional  patterns  which  result  when  n  =  6. 

Variation  in  element  spacing  could  influence  all  of  the 
terms  in  the  equation  for  the  directional  pattern  of  the  array. 
Figures  9  and  10  demonstrate  this  effect,  showing  the  directional 
patterns  which  result  when  the  n^  elements  are  at  spacing  d  , 
while  the  n~  elements  are  at  spacing  2d,. 

In  these  figures  and  in  other  investigations  of  the 
directional  patterns,  it  is  apparent  that  the  cosCdo/d^  )  v|J  term 
caused  by  separation  of  group  centers  has  the  greatest  effect 
in  narrowing  the  principal  lobe. 

C.   Data  Processing  by  Correlation  Circuits. 
Data  processing  by  correlation  circuits  offers  an 
attractive  possibility  for  improving  the  presentation  of  the 

formation  available  from  an  antenna  array.   Such  a  system 
could  be  based  on  the  multiplication  of  the  output  voltage 
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of  one  linear  array  by  the  output  voltage  of  a  second  linear 
array  in  which  a  progressive  phase  shift  has  been  introduced, 
For  example,  using  again  a  total  of  six  elements: 


(UUUUi 

px  0V  $3  <$«  $r 

V     V 3z id b__ 


L 


AV6. 


By  providing  a  number  of  phase  shifted  circuits  with  the  pro- 
gressive phase  shifts  chosen  to  form  directional  patterns 
whose  maxima  occur  at  the  zeros  of  the  unshifted  pattern  or 
at  the  sidelobe  maxima  of  the  unshifted  pattern,  a  number  of 
such  correlation  processes  can  be  carried  out.   By  combining 
these  correlation  voltages  into  a  single  output  voltage,  an 
improved  directional  pattern  is  formed. 

Figure  11  shows  this  resulting  pattern  for  the  six- 
element  array.   The  unshifted  pattern  is  given  by 

Ps(V)    =  sinoC^  12) 
F^Y)        6  sin(V/J) 

The  phase  shifted  patterns  are  given  by 

where  <»*/6,   */3,  ir/2,  2TT/3,  5TT/6,  and  tt  . 
By  fully  utilizing  all  of  the  phase  shifted  circuits  connected 
to  these  six  elements,  it  would  be  possible  to  provide  a  number 
of  antenna  patterns  with  main  beams  displaced  (equally  in  ty   ) 
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throughout  360°,  all  from  a  stationary  array.   This  multiple- 
beam  antenna  is  not  new,  but  the  improved  beamwidth  provided 
by  the  correlators  permits  a  refinement  in  performance  without 
increasing  the  overall  length  of  the  array. 

Using  the  same  six  elements  to  provide  unshifted  as 
well  as  phase  shifted  voltages  will  not  permit  elimination  of 
noise  voltages  by  averaging.   And  it  does  not  take  advantage 
of  the  beam  narrowing  feature  of-  the  cos^do/d}    tern,  present 
when  the  linear  arrays  are  separated.   Thus  some  improvement 
in  performance  would  be  expected  from  splitting  the  elements 
into  two  separated  linear  arrays,  with  progressive  phase 
shifting  applied  to  one  of  these  resulting  arrays.   Figure  12 
shows  the  pattern  which  results  when  the  six  elements  are 
formed  into  two  uniformly  spaced  three-element  linear  arrays. 


I 


I    I    I   I 

4i  I  fer|  lfy 


* 


p   (V\        sin  3(4*72)      sin  3  %(4>t*Q    -      m**^ 
P6(T}   "   5    sin    (V^/ij    3    sin  lfcj±«c  J   cos (3^*) 

where  *C=0,    TT/3,    2TT/3,    and  T*  . 
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APPENDIX  I:   THE  PROBABILITY  DISTRIBUTION  FOR  THE  OUTPUT  NOISE 
VOLTAGE  OF  A  TWO-ELEMENT  CORRELATION  ARRAY. 

\  ran    distribution  of  noise  sources  surrounding  a 
two-element  antenna  array  will  induce  voltages  which  can  be 
represented  by  stationary  time-  series  with  gaussian  probability 
distributions.   The  noise  voltages  car.  thus  be  completely  de- 
scribed statistically  by  specifying  their  cross-  and  auto- 
correlation functions.   Each  noise  source  will  induce  a  voltage 
on  one  element  which  is  identical  to  that  induced  on  the  other 
element,  but  with  a  time  delay  factor  d^«=  to  the  finite  time 
of  travel  of  the  wavefront  frcm  one  element  to  the  other.   As 
shown  previously,  the  cross-  and  auto-correlation  functions 
of  the  total  induced  voltage  can  be  obtained  by  summing  the 
contributions  of  the  individual  sources.   For  a  uniform  distri- 
bution of  independent  noise  sources  and  rectangular  band 
limiting : 

^     r    \        i      o     sinirBs  rl  //>x  sin*irB« 

Tu(s)  =  +  n  :>n  tBs    cos  *,.  s   -  Rn(0)  s^s°"    cos  a3„  s 


qj  f.        D  ...  sin  kSM.  cin  !lM 


sin  irBs 


k  TrBd  irBd  ^£\7    cos<^s 

c  c 


=  r_  cr*  >  sin  it  Bs 

c^K  L  J    COS  CU-S 

ttBs 


where : 

1 )  noise  sources  in  element  of  angle  dft  produce  noise 
power  y%  dXL  watts/cps; 
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2)  B  is  bandwidth,  CO      is  center  frequency; 

3)  k  = —  ,  and  d  is  element  spacing. 

tiB 

The  elements  of  the  system  can  be  represented  as: 

n'(t)  n»(t) 


x 

T 


;w 


■//lien  n,(t)  =  n2(t),  the  circuit  is  that  of  a  square  1 

detector;  the  statistics  of  this  system  have  been  studied 
several  authors  L  *    L*.   Emerson  has  shown  that  the  probability 
distribution  of  the  output  voltage  of  the  square  law  detector 
is  determined  by  the  eigenvalues  of  a  linear,  homogeneous, 
integral  equation  whose  kernel  is  a  function  only  of  the  cor- 
relation function  of  the  input  voltage  and  the  impulse  response 

of  the  circuit  employed  as  an  averager. 

1  ° 
Using  this  general  method  Lampard   has  obtained  an 

explicit  solution  for  correlated  input  voltages,  but  only  for 
the  case  of  no  postmultiplier  filtering. 

Although  it  appears  that  there  is  no  explicit  solution 
to  this  general  antenna  noise  problem  when  postmultiplier 
filtering  is  employed,  the  resulting  distribution  can  be  repre- 
sented by  an  infinite  series  whose  significant  terms  can  be 
calculated.   If  these  input  noise  voltages  have  bandlimited 
gaussian  distributions,  the  output  from  the  averaging  circuit 
can  be  expressed  in  terms  of  an  infinite  series  and  displayed 
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graphically  in  a  manner  similar  to  that  used  by  Emerson. 

The  output  voltage  of  the  circuit  is  written  as  a  con- 
volution 

vQ(t)  =  ^   nL(t-u)n  (t-u)h(u)  du 

where  h(u)  is  the  impulse  response  of  the  averaging  circuit, 
Since,  for  this  antenna  situation,  n,(t)  differs  from  n„(t)  only 
by  the  time  delay  factor  t  ,  the  output  can  be  written 


vo(t)  =  ^   n1(t-u)n1(t+  t  -u)h(u)   du. 


The  moments  of  the  output  distribution  can  be  determined  di- 
recti y : 

vQ(t)  =  \  n1(t-u)n1(t+  X  -u)h(u)  du. 
0 

For  this  perfect  averaging  circuit 

h(u>  At  •  °*utt 

I  0  ,  elsewhere 
Then 


h  «  ^tft  -  \\  Vu  (o) 


du 


where  f  (s)  =  n^(t)n  (t+s)   are  the  correlation  functions  for 
the  bandlimited  noise  voltages. 

K  =  v  (t)z  =  i  C   C  n1Ct-u;n2U-u;n1Ct-v;n2U-v;  du  dv 


l2t     $o[tl2(0>2   +^11(U"V)2    +Vl2<u-v)2]du 


dv. 
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The  variance  is  <T     K0  -  K 


1 


A  general  expression  for  the  moments  could  be  written 
in  the  form  of  an  interated  integral 

K.  =  v^t)1  =  1.  f   .  .  I  nL(t-x  )n2(t-x1)~  -^(t-x^Ct-x^ 

J0      J0  dxL dx.  . 

The  output  probability  distribution  can  now  be  determined  in 
terms  of  a  Gram-Charlier  series,  type  A.   •     If  we  introduce 
the  normalized  variable 

y  =  v°(c)  -  \^    , 

then  a  convenient  expansion  for  the  probability  distribution 
of  y  is  given  by  the  series  grouped  according  to  Edgeworth: 

P(Y)  7ke^y1  1+C3H3(y)  +  Wy)  +  Wy)  + ]  • 

Where  H  (x)  is  the  Hermite  polynomial  of  degree  n  defined  by 


— n  (e-***)  =  (-1)   Hn(x)  e"2X   , 


where 


c,  -  All 

an  t 

and  u  ■    is  the  i-th  central  moment,  obtained  from  K. ,  the  i-th 
moment . 

To  illustrate  the  change  in  the  probability  distribution 
of  the  output  voltage  with  increasing  averaging  time,  let  us 
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consider  two  isotropic  antenna  elements  in  a  uniformly  distri- 
buted noisy  medium  with  an  element  spacing  d  =  ^  .   For  this 

2  "* 

element  spacing,  R  (X  )    =  r   R  (0)« 

n       • 

With  no  postmult  iplier  filtering  the  distribution  is 
given  explicitly  by  La  rv.p  a  re  : 


P(vo)  -  i[Va(o2-Y12(o)^ ' >^° y^oky  12(o>^ 

K4  vu^)2-^^  j  • 


where  K  (x)  is  the  modified  Bessel  function  of  the  second  kind 
o 

and  zero  order. 

With  postmultiplier  averaging  included,  the  equations 
for  the  first  three  central  moments  of  the  distribution  are: 

,   t  o 

/A.  =   v~TtT  =  -(  Rn(X)  du  =  R  (T)  -  i   R  (o). 

'   1     O        tj  n         TTn 

0 
P>2    =<S   =  I  C  (jik(C)2^(X)2lsin27rBW-y?  cosL„(u-v)ldu  dv. 


A 


i  .  _ 


cos  u3#(vu  -  y  )  cos  u»«  (y-w)cos  <*»#  (w  -  u^  )(du  dv  dw  , 

The  output  probability  density  functions  are  now  expressed  in 
terms  of  the  Gram-Charlier  series.  They  have  been  plotted  in 
Figure  13  for  selected  values  of  the  bandwidth- averaging  time 
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product,  Bt .   The  coefficients  have  been  obtained  by  neglecting 
the  terras  in  the  cosuJ  (u-v)  products  which  represent  the  high 
frequency  residue.   Since  in  a  typical  circuit ^ 0>>   B,  they 
are  vanishingly  small  for  values  of  Bt  considered  here.   Actual 
solution  of  the  integrals  to  obtain  the  series  coefficients 
was  carried  out  on  a  Bendix  G-15  digital  computer. 

For  the  case  Bt  =  0  (no  averaging)  the  curve  is  exact, 
and  the  effect  of  all  frequencies  is  considered. 

For  Bt  =  1 ,  Bt  =  2 ,  and  Bt  =  4,  the  convergence  is  fairly 
rapid  (three  terras  of  the  series  required)  for  graphical  ac- 
curacy in  regions  near  the  mean.   Accuracy  for  the  extremes  is 
not  claimed,  but  would  require  additional  terms  in  the  series. 

It  will  be  noted  that  the  tendency  toward  normality 
with  increasing  Bt  (increasing  average  time)  is  pronounced. 
This  is  in  general  agreement  with  the  results  observed  by 
Emerson  for  the  square  law  detector  and  would  be  expected  from 
the  general  behavior  of  integrated  noise  derivable  from  the 
Central  Limit  Theorem. 


73, 


'   fi 

<D 

rf 

B 

m 

0) 

t— 4 

.— < 

<u 

o 

* 

00 

rt 

00 

*-> 

u 

"o 

E 

> 

4-1 

0) 

(A 

M 

•r-> 

c 

o 

■  r-1 

C! 

■t-> 

3 

4) 

«*' 

3 

> 
it 

^ 

0 

0) 

II 

OJ 

S 

XI 

XI 

C 

<-> 

-  r~* 

W) 

Ih 

c 

O 

> 

u 

V-i 

^ 

nJ 

10 

O 

a 

C 

<4H 

CO 

o 

■tH 

>- 

■»-> 

4-> 

(4 

c 

u 

*H 

0) 

C 

"3 

6 

>> 

0 

I— 1 

W 

•iH 

H 

(A 

•»-> 

C 

n3 

0) 

t; 

"0 

>~ 

>- 

b 

+-> 

c 

•rt 

u 

•— i 

■r* 

XI 

aj 

J3 

O 

U 

i 
i 

o 


in 

d 

i 


74. 


APPENDIX  II: 


CALCULATION  OF  VOLTAGE  RESPONSE  PATTERNS  WITH 
TWO  SIGNAL  SOURCES. 


The  antenna  will  be  a  four-element  correlation  array  with 
relative  spacings  d,~  =  D,  (I23  =  D,  cL,  =  2D,  and  D  =  ^/2. 
The  voltage  induced 


Ugh^ 


Z 


*x 


*t 


on  the  i-th  element  by  signal  source  A  is  V  =  A  coscu   (t  +X       ) 

A  A       Ai 

and  that  by  signal  source  B  is  VR  =  B  cosws(t  + 1      )  y    where 
^Ai  =  (di/c)sin  0A  and  tT   =  (di/c)sin  0„.   The  signals  which 

are  interesting  in  this  study  of  correlation  processes  are 

grouped  into  two  cases:  (1)  coherent  signals  (**^A  =  ***n)  and 

A      B 

(2)  signals  with  slightly  different  frequencies  (U>  „-  u)  =  $u>). 

Therefore,  (u)  d-/c)*(ci)  d./c)  and  the  following  notation 

A  *-  D   i- 

can  be  used : 


X  -L£  sin  *A 


and 


y  =*£_  sin  0 

c       B 


Then  the  product  of  the  four- element  voltages  can  be  written: 


I' 


=\  [A  cos  cx>At   +   B   cos"JRt][  A  cos    ( o>At   +  X)    +   B   cos    (^3t    +  Y)]. 


[A  cos(toAt   +   2X)   +   B   cos    (tOBt    +   2Y)][A  cos    (cuAt+4X)  +  B  cos  («ty4 

Carrying  out   the    indicated  multiplications,    and   rejecting  the 

< 


i 
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multiple  frequency  components  (2  oj   and  4u>),  ^ives 

V  =  A^[cos  X+  cos  3X+  cos  5X]  +  ^[cos  Y+  cos  3Y+  cos  5Y] 
+  A^B[Cos(x+S«t)+cos(3X-J«t  )+cos(5X+S«t)+cos(X+4Y+ln»t) 

+  cos(X-4Y-Sa>t)+cos(2X+Y+S«t)+cos(2X-Y-Stt>t)+cos(3X+2Y+S*t) 
+  cos(3X-2Y^t)+cos(3X-4Y-Su>t)+cos(5X-2Y-&-t)+cos(6X-Y-&ut)  ] 
+  i^[cos(Y-Su>t)-t-cos(3Y+<k>t)+cos(5Y+£uit)+cos(Y+4X-S*t) 

o 

+  cos(Y-4X+£u*t)+cos£Y+X-i*t)+cos(2Y-X+Ju>t)+cos(3Y+2x4tt/t) 

+  cos(3Y-2X+iuit)+cos(3Y-4X+W)+cos(5Y-2X+Suit)+cos(6Y-X+^t)] 

2  2 
+  ^-[cos(X+2Y)+cos(X-2Y)+cos(X+4Y)+cos(X-4Y)+cos(2X+Y)  +  co^2C-Y) 


+  cos(2X+3Y)+cos(2X-3Y)+cos(3X+2Y)+cos(3X-2Y)+cos(4X+Y) 
+  cos(4X-Y)+cos(X-6Y-2iut)+cos(2X-5Y-2iVjt)+cos(3X-4Y-25iut) 
+  cos(4X-3Y-2S"wt)+cos(5X-2Y-2&>t)+cos(6X-Y-2&t)] 

Although  this  is  a  formidable  expression,  the  capability 
of  the  array  to  resolve  the  two  sources  can  be  estimated  by 
examining  two  simplified  equations: 
Case  1.  .  .when  source  A  is  located  on  the  principal  lobe  of 

the  array  (X  =  0), 
Case  2.  .  .when  the  principal  lobe  of  the  array  is  directed 

between  the  two  sources  (X  =  -Y). 
If  it  is  assumed  that  A  =  B  and  that  «u>  =  0,  then  the  equations 
for  these  two  cases  are: 

Case  1 

a4 
V  =   ^.[6   +    12    cos   Z'    ♦   10    cos   2Z'    +  8    cos    3Z1    +   6   cos    4Z' 

+  4   cos   5Z'    +   2    cos   6Z' J 
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Case  2 


V  =  £L-[12  cos  Z"  +  12  cos  3Z  "  +  12  cos  5ZM  +  12  cos  7ZM  j 
E 

where,  with  the  angular  displacement  of  the  sources  denoted  by  0, 

Z'  =  TT  sin  0,     and     Z"  =  tt  sin® 

Resolution  of  the  two  sources  is  certainly  possible  when  the 

magnitude  of  the  voltage  response  with  the  main  lobe  directed 

between  the  sources  (case  2)  is  less  than  the  magnitude  of  the 

voltage  response  with  the  main  lobe  directed  at  one  of  the 

sources  (case  1).   Although  this  is  not  the  only  condition  under 

which  resolution  will  occur,  a  comparison  of  these  two  voltage 

responses  as  functions  of  the  separation  of  the  sources  should 

indicate  the  general  action  of  the  array  in  the  presence  of  two 

sources.   Resolution  is  certain  for  separations  of  approximately 
1.0 


Angular  separation  of  sources,  0,  in  degrees 

VOLTAGE  RESPONSE  OF  ARRAY  AS  A  FUNCTION  OF  SEPARATION 
OF  SOURCES 

Dashed  line:   main  lobe  directed  at  one  source 
Solid  line:   main  lobe  directed  between  sources 


77 


14  ,  29  ,  and  43°  where  the  voltage  response  for  the  main 
lobe  directed  between  the  sources  drops  to  zero.   It  is  interes- 
ting to  note  that  the  voltage  response  when  the  main  love  is 
directed  at  one  of  the  sources  is  uniformly  low,  and  the  array 
will  not  indicate  a  maximum  voltage  response  for  this  orien- 
tation.  The  general  voltage  response  patterns  of  this  array 
for  several  specific  source  separations  are  shown  in  the 
accompanying  diagrams,  (Figs.  14,  15,  and  16). 
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FIG.    14.    --  Voltage  response  pattern,    four-element  correlation 
array  and  two  coherent  signal  sources.      Voltage  normalized  to 
unity  for  maximum  response  of  two  coincident  signal  sources. 
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Separation  of  Sources  =  34 


Source  B 


FIG.  15  --  Voltage  response  patter,  four-element  correlation 
array  and  two  coherent  signal  sources.  Voltage  normalized  to 
unity  for  maximum  response  of  two  coincident  signal  sources  . 
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FIG.    16  --   Voltage  response  pattern,    four-element  correlation  array 
and  two  coherent  signal  sources.      Voltage  normalized  to  unity  for 
maximum  response  of  two  coincident  signal  sources. 

A  similar  response  pattern,    resolving  the  two  sources,    can  be 
obtained  with  a  linear  antenna  with  aperture  3X.    (6D).      This  is 
approximately  the  same  length  as  the  correlation  array  whose 
patterns  are  plotted  here. 
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APPENDIX  III:   APPLICATION  OF  CORRELATION  TECHNIQUES  IN 

LINEAR  ARRAYS. 


A.   INTRODUCTION 

The  consideration  of  che  statistical  nature  of  noise 
and  the  application  of  correlation  techniques  can  be  utilized 
in  the  study  of  antenna  systems  which  are  designed  to  receive 
a  signal  embedded  in  a  general  noisy  environment.   The  tech- 
niques of  correlation  theory  can  be  employed  to  indicate  the 
method  of  operation  o£  the  antenna  system  which  will  optimize, 
in  some  specified  sense,  the  reception  of  the  desired  signal. 
It  is  the  purpose  of  this  report  to  summarize  some  previous 
results  in  this  field  and  to  indicate  a  general  methodology  for 
attacking  the  problem. 

The  general  case  considered  is  that  of  a  single  signal 
source,  a  point  source,  in  a  noise  medium.   No  specific  restric- 
tions are  imposed  on  the  type  of  signal  considered;  only  a 
statement  of  total  signal  power  is  made.   Thus,  the  treatment 
is  generally  applicable  to  a  signal  which  is  a  random  function 
having  the  same  spectrum  as  the  background  noise.   A  signal  which 
is  correlated,  or  whose  spectrum  varies  from  that  of  the  back- 
ground noise  would  obviously  permit  further  modifications  of 
the  receiving  techniques  which  could  give  improved  antenna 
system  performance. 

The  receiving  apparatus  considered  is  a  linear  broad- 
side array,  oriented  to  put  the  signal  source  at  the  point  of 
maximum  sensitivity.   Under  these  conditions  the  signal-to-noise 
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ratio  of  the  combined  output  ot  the  array  car  be      nized 
with  respect  to  both  of  these  variables. 

The  application  of  correlation  techniques  to  receiving 
systems  has  been  considered  in  many  forms  in  the  literature. 
The  present  method  of  analysis  is  based  on  that  introduced  by 
Faran  and  Hills^  in  their  studies  of  acoustic  systems. 


Each  element  of  an  array  of  antennas  in  a  noise  fii 
can  be  considered  as  an  individual  source  whose  output  can  be 

-tably  modified  in  amplitude,  phase,  or  time  delay  before 
beinu  combined  with  other  source  outputs  to  prouuce  the  desired 

;al  effect.   The  output  of  each  such  element  will  be  a  de- 
sired signal  voltage  plus  an  undesired  noise  voltage .   If  there 
is  a  correlation  between  the  noise  voltages  occurring  on  any 
two  of  the  antennas,  then  it  should  be  possible  to  modify  and 
combine  these  voltages  in  such  a  manner  as  to  reduce  the  total 
effect  of  the  noise  in  the  output:  of  the  system. 

The  correlation  coefficient  used  will  be  the  time  average 
function  defined  by: 

RL2(T)  =    lim   ^A       ei(t)e2(t+T)dt 

-T 


When  ei(t)  and  eo(t)  represent  the  instantaneous  voltages  on 

2     — ? 
two  antennas,  e]_   and  e.^      are  proportional  to  the  average 

powers  of  the  two  sources.   When  these  two  voltages  are  added, 
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the  average  power  of  the  comgination  is  proportional  to: 

e  =  [ei(tl  )  +  e2  (t2)  I2 


7 


=  eL(t  )fc  +  e2(t2r  +  :.e1(t1)e2(t2) 


(2.2) 


=  e1(t1)-  +  e2(t0)^  +  ^R12(tl  _  tl') 


=  eL(tL)^  +  e2(t  )'   +  2R12(T) 

The  cross-correlation  coefficient,  R   (T)  ,  will,  in  general, 

take  on  both  positive  and  negative  values  as  XT  is  varied. 

This  argument  can  be  extended  to  the  case  of  n  sources. 

In  this  case  the  power  of  the  combined  output  is  proportional 

to: 

n  n     n 

7=Z  ei(v2  =  Z  Zl  e^w 
i=i  i«i  j=i 

which  will  contain  [as  R^CT!)]  the  correlation  coefficients 
between  each  antenna  and  all  others. 

Consider  an  isotropic  antenna  located  at  the  origin  of 
the  spherical  coordinate  system,  with  noise  sources  distributed 
on  the  surface  of  the  enclosing  sphere.   Then  the  noise  arriving 
from  a  differential  angle  d<\   in  the  direction  oC  contributes 
an  amount  f(oc)dot    to  the  average  noise  power  at  the  output 
of  the  antenna.   It  will  be  assumed  throughout  this  report 
that  noise  signals  arriving  from  different  directions  are 
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-  idependent.     e  spectral  density  of  the  noise  power  is  then 
expressed  as  a  function  of  position  and  frequency: 

P  (0,  ^,oO  )  =  I  (0,  4    )  •  W  (a*)  (2.4) 

The  actual  antenna  will  not  be  isotropic  but  will  have  a  charac- 
teristic power  gain  pattern,  G  (0,  (f>) ,  and  the  spectral  density 
of  the  noise  power  [W(t*>  )  ]  will  be  modified  [to  .W  (•»)  ]  by  the 
characteristics  of  the  antenna  and  its  associated  circuitry 
(amplifier  and  phase  shifting  network). 

Since  the  correlation  function  (autocorrelation  function 
in  this  case  of  a  single  antenna)  and  the  power  spectral  den- 
sity are  Fourier  transform  pairs,  the  process  of  determining 
this  correlation  function  is  that  of  the  triple  integration: 

G(9,<rf)W  («)eJ   sin  OdOd^du  (2.5) 


This  method  of  calculating  the  correlation  coefficient 
t  readily  modified  to  apply  to  the  case  of  an  array  of  in 
antennas.   If  it  is  assumed  that  each  of  the  antennas  in  the 
array  is  identical  to  each  of  the  others,  and,  further,  if  it 
is  assumed  that  no  mutual  coupling  exists  between  antennas, 
then  it  is  only  necessary  to  add: 

(1)  a  term  to  account  for  the  phase  difference  (or  time 
delay)  in  the  voltages  induced  by  any  given  noise  source  be- 
cause of  the  distances  between  antennas,  and 

(2)  constants  A^  and  A.  which  describe  the  magnitude  of 
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the  gains  introduced  by  the  amplifier  in  each  antenna  circuit. 

A.. 


ttjct)- VjR,u«>--ySS  1"" 


R, 

J      i  J  J- J  ""   -^r  JJJ 

(2.6) 

Let  the  noise  voltage  received  by  the  i-th  antenna  be 
denoted  by  n.(t).   Then  the  total  noise  power  of  the  array  is 
proportional  to: 


rr. 


NrX  Aini  (t-Ti)" 


(2.7) 
i=l 


where  A.  is  the  gain  introduced  at  the  i-th  antenna,  and  t-    is 
the  time  delay  at  this  antenna. 

mm  mm 

1-1  j  =  l  1=1  j=l 

In  order  to  facilitate  interpretation  of  the  mathematical 
solution  of  the  problem,  it  i.s  convenient  to  normalize  power 
equations.   Since  the  antenna  elements  have  been  assumed  to  be 
identical,  R1 •  (T)  depends  only  on  the  time  delay  T,  and  R1 (0) 
is  the  same  for  all  elements.   Since 

R'(0)  =  lim  ^  \  n2(t)dt  =  nTt72 
T->  oo     {T 

this  is  a  measure  of  the  power  induced  on  a  single  antenna. 
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Then,    for   the    syste 

m        m  R«       (X     -X) 

^~  ^r —  11      i        i 

= y>      y>       AjAj  R'(0)  (2.9) 

1=1      j=l 

For  a  Linear  broadside  array  with  a  single  signal  scarce 
Located  at  the  point  of  maximum  sensitivity  of  the  array,  there 
is  no  time  delay  in  the  reception  of  the  signal  voltages  by  all 
elements, so  signal  power  is  proportional  to: 

m    m 
Sn2  =^>   ^>   AtAj  =  (AL  +  A2  +  .  .  .  +  Ani)2       (2.  10) 
1=1   j=l 

If,  furthermore,  An  +  A^  +  .  .  .  +  A  =1,  then  Eq.  (2.9) 

12  m  n   *    ' 

o 

for  N    is  the  inverse  of  the  signal- to-noise  ratio  of  the  com- 

n 

^d  output  of  the  array. 


C.     TMUM  PERFORMANCE 

For  the  purpose  of  this  report,  optimum  performance  will 
be  that  which  gives  a  maximum  signal-to-noise  power  ratio  in 
the  output  of  the  system.   Variables  to  be  considered  are: 
(1)  the  relative  gains  of  the  amplifiers  connected  to  each 
antenna  (A^),  (2)  a  time  delay  in  each  antenna  circuit  (T.  ) ,  and 
(3)  antenna  spacing  (d). 
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First,  it  is  necessary  to  examine  the  nature  of  the 
correlation  coefficient,  R-  •  (T)  .   Consider  a  linear  array  of 
m  identical,  evenly- spaced  isotropic  antennas. 


Assume  a  continuous,  uniform  distribution  of  white  noise  sources 
over  the  surface  of  an  enclosing  sphere. 
Then 

G(9,  tf)  =  L/41T 

I  (9,  <0  =  1/4IT 

W(«o  )  =  \   watts/cps 

2 

Let  each  amplifier  be  band  limited  vrith  identical  rectangular 

*ands,  UJ   +  &6l>  . 
0 


-u> 


C    2lT   w 


R'^CT) 


64^ 


-OJ   -4*1  0   ''C 


-J— —  cos  0  iw-v 

sin  0  e   c       eJ 


dQdtfdco 


u)  +  &<* 


V*  =JL 


27t  ^tr      «d  cos  e 

sine  e  eJ    dQ  d*  d^ 


0  ^0 


(3.1) 


~(jj 


,U)0  +  4W 


sin«ud/c  +  e-^^  ^  +\   sinwd/c   K   ((3.2) 


w  -fcto  «*d/c 
0 


Ou 


o> 


d/c 
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This    is    solved    to   give 


-  s 


+  s. 


'CO  n  «-AO^ 


*>n 


ZTTd    -a)  1)1  Sj 


Under  the  condition  Acu«fr>^,  this  simplifies  to 


(3.3) 


R'  (d>   T?   j^    sin    2-lTd/fr     cosCO  X 
R'(0)  2nd/A  0 


(3.4) 


Thus,  the  correlation  coefficient  can  take  on  both  positive 
and  negative  values;  then,  if  the  magnitude  of  the  correlation 
coefficient  is  established  by  the  element  spacing,  a  time  delay 
network  can  be  utilized  to  obtain  an  effective  correlation  co- 
efficient  of  any  lesser  absolute  value.   In  this  case,  of  course, 
the  effect  of  this  time  delay  on  the  desired  signal  voltage 
must  be  considered  in  evaluating  the  system  performance. 

Maximizing  the  output  signal-to-noise  ratio  with  respect 
to  the  gain  introduced  at  each  antenna  amplifier  or  with  respect 
to  element  spacing  is  a  matter  of  applying  the  Lagrange  mul- 
tiplier method  to  Eq.  (2.9)  subject  to  appropriate  constraints 
on  total  amplification  or  total  array  length.   The  desired 
quantities  are  found  by  solving  a  set  of  m+1  simultaneous  linear 
equations  in  m+1  unknowns  (for  an  array  of  m  elements).   Sym- 
metry can  be  used  to  reduce  the  number  of  equations  and  unknowns 
considered. 
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D.   EXAMPLE:   FIVE- E      C  ARRAY;  ISOTROPIC  ANTENNAS 

Assumptions : 

1)  Five  identical  isotropic  antennas  in  a  linear  broad- 
side array. 

2)  Single  point  signal  source  located  in  the  direction 
of  the  maximumnof  the  array  radiation  pattern. 

3)  White  noise  sources  uniformly  distributed  over  the 
surface  of  an  enclosing  sphere. 

'+)   No  mutual  coupling  between  antennas. 
5)   Antenna  amplifiers  have  identical  rectangular  band 
characteristics  uj  +6u>,  and  A"**  =  0.01  <**  n . 


V        V        V        V 


A3 


A< 


T, 


OUTPUT 


Case  I. --Uniform  spacing,  (d),  between  antennas. 


4.\L2  +  A22  +  A32  +  A,2  +  A52  +  2R10(d10)A1A, 


12s  12'  12 


+  2A1A3R13(dI3)  +  2AlA4Ru  (d^)  +  2A1A5Ri5(di5) 

+  2A2A0R23(d?3)  +  2A2A4R2,(c:o;)  +  2A2A5R25(d25) 
+  2A3A4R34(d34)  +  2A3A5R35(d35)  +  2A4A5R45(d45) J 


(4.1) 
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and 

2 


Sn  (Ax    +   A,    +   A3    +   A4   +   A^)2    -    1 


rnr   a   symmetrical   array    (A  =A    ,    A  =A   ),    and  with  Ti^T^T^l^^r^O 

Nn2    =   A32    =    2A22    +    2AL2    +    4A1A2R1?(d)    +   4A2A3R23(d) 

+   4A1A3Rn(2d)    +   4A1A?Rl9(3d)    +   2A22r      (2d) 

+■  2A12Ru(4d) 

range  multipliers  give  four  equations  in  four  unknowns: 


>     1      + 

Ssn2 

=   0 

w  + 

S  A2 

*  A2 

=   0 

o  Nn2      ^  V 


+   — —  =  o 


°    3         3 


S  2  =  1 

n 


The  correlation  coefficient  as  a  function  of  element 
spacing,  d,  can  be  developed  as  in  Section  III. 
G  (O  ,<p   )  =  l/4lT 

I  (&,*  )  =  1/47T 

W(u>  )  =  ^-watts/cps 
T  =  0 
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'(0,    0)         ?TT       AJ"1  ~ 


(  CO  n+  Ac*) 


■  type?  j"..  r  L.0]  M0d  I 


Again,    in    some    c    •        ,    this     sine-inti      ral    expression   can   be 
tec     by 

s  in  co  <  ' ; 

l 

^  c«':.'l- 
The  two  expression      equal  for  .it  Least  twe     nific< 
This  is  show     Fig.  ' -7 • 
For  s       (d)  less  than  a  ;th,  the  a 

fications  of  indivic1     mtenna  voltages  which  will  j     t  in 
a         signal-to-noise  ratio  in  tue  combined  output  differ 
jnificantly  from  the  uniform  amplification  system.   For  . 
s  of  a  half-     e  ter,  optimum  amplification 

tt       '.     resu        the    ial         »ve  give  only  a 
minor  improvement  in  signal-to-noise  ratio  over  that  obtaine 
with  ..r.iform  amplification.   Figure  2  shows  this  signal-to- 
noise  ratio  as  a  function  of  element  spacing  for  the  optimu 
lification  case  and  for  the  uniform  amplification  case. 
Two  observations  can  be  made  from  an  examination  of 
Fig.  IS: 

1)   For  this  case  of  uniform  spacing,  maximum  signal- 
to  noise  ratio  occurs  for  a  spacing,  d,  of  315   between  elements 
Reference  to  Fig.  1  shows  thai  for  this  spacing  the  cross- corre- 
lation coefficients  between  all  antennas  (d,  =  315°,  d   =630°, 
d->  =  945°,  d,  =  1260°)  are  either  negative  or  negligibly  small. 
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Optimum 
amplification 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


Uniform 
amplification 


+ 


4- 


+ 


90' 


18(T  270u  360° 

Element  spacing  (in  degrees) 


450 


FIG.  18  --  Five-element  array,  uniform  spacing,  isotropic  antennas, 
single  signal  source,  spherically  distributed  noise  (0  db  corresponds 
to  signal=-to-noise  ratio  of  unity  on  a  single  element). 
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,  t  2.   effect  of  the  cr  ents  of  noise  power  is 

to  decrease  the  noise  power  Inter    t ed  by  each  of  the  five 
antennas  individual!;  . 

2)   A  significant  incre    In  signal-to-noise  ratio  can 
be  attained  at  close  spacings  '     ilizing  the  amplification 
patterns  resulting  from  a  solution  of  Eq,  (4.4).   The  rapid 
increase  in  individual  antenna  voltage  amplifications  for  op- 
timum performance  as  the  spacing  is  narrowed  is  shown  in  Fig. 
and  listed  in  detail  in  Table  L.   Calculations  for  spacing 
in  this  region  involve  small  differences  of  large  quantities 


Element 
Spacing 

A 

if ications 

$L 

A  2 

A3 

A4 

i_5 

30° 

+  57.5 

-212. 

311.7 

-212.9 

57.5 

45° 

13.3 

-  44.9 

64.2 

-  44.9 

13.3 

,   50° 

4.03 

-11.56 

16.1 

-11.56 

4.03 

75° 

1.81 

-4.197 

5.77 

-4.197 

1.81 

90° 

.98 

-  1.68 

2.40 

-  1.68 

0.98 

12C° 

0.417 

-0.214 

0.595 

-0.214 

0.417 

150° 

0.253 

+0.116 

0.264 

+0.116 

0.253 

180° 

0.199 

+0.201 

0.199 

+0.201 

0.199 

TABLE  1. --Individual  antenna  voltage  amplification  to 
produce  maximum  signal-to-noise  rationfor  the  given 
element  spacing. 

so  the  sine  integrals  cannot  be  approximated  by  the  s^-n  x 

function.   As  the  spacing  becomes  small,  optimum  performance 

depends  on  precise  adjustment  of  the  gains  of  the  individual 
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Element  spacing  (in  degrees) 


FIG.    19  --   Five-element  array,    uniform   spacing,    amplification  of 
antenna  voltages  to  produce  maximum  signal-to-noise  ratio. 
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lplifiers,    and,    there  increasingly   dif- 

:uLt    to    att  from   a   practic  nt   of   view.      This    is 

Logous    to    the    problei  ing    in    super-gain 

antennas. 


^7 

— y — 


V 


& 


Ax 


w 


?3 


As 


h  se  II.-  Variable  spacing  between  antennas. 
Co:  Isotropic  elements  in  a  symmetrical  linear 

brc     de   array  with  unequal  ele  ent  spacing.   Assuming,  as 

that    =  A5,  A2  =  A4>  and  XL=X?=T3=X4=t5=0 ,  the 
equations  del      noise  power  and  signal  power  become: 
Nn2  ={A3        ^V  +  4A2A3R23(x)  +  4A1A2R12(y) 

+  '       j(x+y)  +  4A1A2R12(2x+y)  +  2A-,  2RU  (2x+2y) 

V   .   /)  (4.6) 


n 


*  -.         J 
=  (A_  + 


ct\~)       -    i 


(4.7) 


it  is  furt       umed  that  overall  length  of  the  array 
is  great  en     to  permit  approximation  of  the  sine  integral 

.. actions  by  — ind  Co  permit  use  of  a  uniform  gain  pattern 

on  the  antenna  amplifiers,  the  Lagrange  multiplier  method  can 
be  readily  err       to  obtain  optimum  spacings. 
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A]_    =    A.    =    A3    =    A4   =    A5    =       .  (4#8) 

=   0.20   +  0.16sin   x   +   n   ifisin  y  +  0tl6sin    (x+y) 
n  x  y  (x  +  y) 


y 

(4.9) 

+  0.16sii](?x'fY)  +  0.08sif)(2x)  +  n  nP.sin(2x+2y) 
zx+y  zx  Zx+Zy 

The  constraint  that  must  now  be  imposed  is  the  total  leagt1 
the  array: 

x  +  y  =  K  (4.10) 

These  conditions  lead  to  three  equations  in  three  unknowns: 


£  ••  2 

O     n 
Sr 

+ 

A 

TIT 

r 

ov 

+ 

* 

5  k 

T7 

=  C  (4.11) 

w  y 

x  +  y  =  K 

These  are  readily  solved  to  give  the  following  criteron  for 
extrema  of  the  signal-to-noise  ratio  as  a  function  of  element 
spacing: 


=  0  (4.12) 


(sin  x\   -   /sin  y\    +  /  ein(2x=y)A   +   /sin  2x1 _ 
x  yl     ^  y  J  ^   2x+y  )  \^J~ 

where  the  prime  indicates  the  first  derivative  of  the  function 
with  respect  to  the  total  argument.   In  general,  roots  of  this 
at  ion  will  give  spacings  for  minima  as  well  as  maxima  of  the 
Lgnal-to-noise  ratio,  so  it  is  necessary  to  carefully  select 
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the  correct  roots.   This  equation  is,  of  course,  a  specific 
application  of  a  general  equation  governing  a  five-element 
linear  array  which  can  be  expressed  as 

R' (x)  -  R'(y)  +  R'(2x+y)  +  R» (2x)  ■  0 
where  R' (  )  is  the  derivative  of  the  normalized  correlation 
coefficient  with  respect  to  its  total  parameter. 

The  maximum  signal-to-noise  ratio  obtainable  by  varying 
element  spacing  as  a  function  of  the  total  array  length  is 
shown  in  Fig.  20.   Spacings  which  give  this  maximum  ratio  are 
shown  for  several  array  lengths  in  Table  2.   Again  it  is  evi- 
dent that  optimum  spacings  results  from  cross-correlation  co- 
efficients which  act  to  reduce  the  total  noice  power  of  the 
array. 


Array  Length 

Spacings 

X 

y 

Signal-to-noise 
ratio  (db) 

360° 

172.0° 

188.0° 

7.086 

450° 

215.5° 

234.5° 

7.891 

540° 

285.5° 

234.5° 

8.511 

630° 

320.0° 

310.0°^ 

8.863 

720° 

282.5° 
492.5° 

437.5° 
227.5° 

7.536 
7.320 

310° 

244.0° 
562.0° 

566.0° 
248.0° 

7.594 
7.759 

900° 

— 

281.5° 
646.5° 

618.5° 
253.5° 

8.025 
8.083 

TABLE  2. --Element  spacings  giving  maximum  signal-to-noise 
ratio  for  various  total  array  lengths  (note  two  possible 
solutions  for  spacings  of  720°,  310°,  and  900°). 
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FIG.    20  -  -   Five-element  array,    variable   spacing  for  fixed  length, 
spacing  adjusted  for  maximum   signal-to-noise   ratio. 
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The  curve  passes  through  a  definite  maximum  which  corresponds 
in  total  array  length  and  in  resulting  signal-to-noise  ratio  to 
the  optimum  condition  discussed  in  Case  I  and  shown  in  Fig.  2. 

In  both  cases  considered  in  this  example,  solution  of 
equations  developed  in  Sections  II  and  III  under  the  conditions 
of  the  assumed  model  have  led  to  definite  conditions  on  the 
array  length,  element  spacing,  or  antenna  voltage  amplification 
which  result  in  a  maximum  of  the  signal-to-noise  ratio  in  the 
output  of  the  total  system. 

E.   EXAMPLE:   FIVE  -ELEMENT  ARRAY;  SHORT  DIPOLES 


Retain  the  same  assumptions  made  in  the  previous  example 
with  a  change  from  isotropic  antennas  to  short  dipoles.   The 
dipoles  can  be  oriented  with  their  axes  perpendicular  to  the 
line  of  the  array  or  parallel  to  the  line  of  array.   Both 
cases  will  be  considered. 


1(0,4*  )  -  l/4r 

G(G,f  )  ■  |(1  -  sin2gCos2t  ) 

W(co  )  *  i- 


Case  I. --Axes  of  dipoles  at  right  angles  to  the  line 
of  the  array. 


R'.     (d,X)  =     l 


n 


0 
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jjeiUosQ   ^ 

3/  ?  ?    x        c  i*/l* 

2(l-sin"0cos^^)e  e1      sLnOdOd^d. 


U»        +Att» 


2IT< 


(5.1) 


2(l-sin20cos2^)e 


3^±   cos© 


^ 


sinO   dO  d^dc*^ 


u) 


0  0      0 


5 


Normalizing  and  integrating  gives: 


R'.  ,(d,t) 
R»(0) 


0 


sin. 


44)    d 

0 


3/v3/2\^— y 


(cos  «^0X) 


(5.2) 


w.ere  ^  3/2i^~) 


is  the  Bessel  function  defined  by 


A  JU  (X  > 

-A-V(x)  =  r(v»+D 


(!) 


In  this  integration  it  has  been  assumed  that  the  rectangular 
pass  band  of  the  amplifiers  is  so  narrow  that 


cu  +AW 


f  (o»  )  dou  «■  f  (tu  )  AtO  . 

0 


*Jr 


The  normalized  correlation  coefficient  has  been  plotted 
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for  T=0  and  for  variable  spacing,  d,  in  Fig.  21. 

The  signal-to-noise  ration  (in  db)  of  this  array  of 
short  dipoles  is  plotted  as  a  function  of  element  spacing  in 
Fig.  22.   As  spacing  is  increased,  this  signal-to-noise  ratio 
passes  through  a  maximum  corresponding  again  to  the  spacing 
which  results  in  a  negative  correlation  coefficient  between  the 
noise  voltage  on  any  antenna  and  the  noise  voltages  on  all 
other  antennas  of  the  array. 


I  (0,  4)   =   l/4ir 

G  (0,  *0  =  1   sin2  0 


W  («o  )  =  A 

2 


Case  II. --Axes  of  dipoles  in  line  with  the  array. 

2i  jr 


*      1     V  \        I     3      .    3.    J!¥  cose 


R'.^d.t)^     \  \       l|sin3ee      c     ^ej*tded^du, 


(5.3) 


<*>      +4o) 


J  ill 

0  0      0 


-J—   cos   0        t 
sin3©   e  eJ    *  d©dgfd"> 


Normalizing  and   integrating,    again   taking  advantage  of  tm  narrow 
bandwidth   to   approximate   the   integration  with  respect   tout,    gives: 

R'^Cd.T)  A  /^0d\ 

— i3 =    JL  J     coswQX  (5.4) 

R'(0)  3/2^    c   J  ° 
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FIG     22  --Five-element  array,    uniform  spacing,    signal«tO- noise  ratio 
as  a  function  of  element  spacing;   short  dipoles  with  axes  P«P«n^cular 
to  line  of  array;  0  db  corresponds  to  signal- to-noise  ratio  on  a  single 
antenna. 


540 


105. 
where,  again 


J&- ~- 


_A_  (^oi\=  r(5/2)  21A 


m 


This  normalized  correlation  coefficient  is  shown  in 
Fig.  23.   The  signal-to-noise  ratio  (in  db)  of  this  array  of 
short  dipoles  is  plotted  as  a  function  of  element  spacing  in 
Fig.  24.   In  this  case,  maximum  of  the  signal-to-noise  ratio 
is  not  pronounced  and  is  only  slightly  greater  than  that  which 
occurs  when  the  noise  is  completely  uncorrelated.   This  could 
be  anticipated  by  observing  in  Fig.  23  the  rapid  decrease  in 
magnitude  of  the  correlation  coefficient  with  element  spacing. 
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FIG.    24  --   Five-element  array,    uniform   spacing,    signal-to- noise  ratio 
as  a  function  of  element  spacing;   short  dipoles  with  axes  parallel  to  line 
of  array;  0  db  corresponds  to  signal-to-noise  ratio  on  a  single  antenna. 
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